$ M1004 : THE-OPEN UNIVERSITY e 
E A : 

ar ¿ : 

a DE Curso Básico de Matemáticas Unidad 4 


Diferencias finitas 


The Open University 


Curso Básico de Matemáticas Unidad 4 


DIFERENCIAS FINITAS 


Preparado por el Equipo del Curso Básico de Matemáticas 


Traductor: 
Hernando Alfonso, Profesor de Matemáticas 
Universidad Pedagógica Nacional 


libros mcgraw-hill 


méxico panamá madrid bogotá sao paulo nueva vork 
londres toronto sidney johannesburg 


dusseldorf singapur 


Queda terminantemente prohibido reproducir este libro, total 
o parcialmente, sin autorización escrita del editor 


RESERVADOS TODOS LOS DERECHOS (D. R.) 


Copyright O 1974 respecto de la edición en idioma español, por 
EDITORIAL McGRAW-HILL LATINOAMERICANA, S. A. 
Colón, República de Panamá 


0-07-090764-1 


Traducido de FINITE DIFFERENCES 
Copyright O 1971 by the Open University Press 
Bletchley, Bucks, England 


IMPRESO EN COLOMBIA PRINTED IN COLOMBIA 


El concepto de enseñanza abierta desarrollado por la Open University 
utiliza materiales en forma de libros, películas y cintas magnetofóni- 
cas que el estudiante utiliza fuera de la institución. Sin embargo, es- 
tos materiales también pueden utilizarse dentro de programas tradi- 
cionales. Los componentes de los programas desarrollados por la Open 
University son: 


LIBROS: Todos los libros tienen el mismo formato, y fueron escri- 
tos, diseñados y producidos por especialistas de la Open University. 
La mayor parte de los libros están profusamente ilustrados con dia- 
gramas y fotografías —muchas a colores— y tienen un diseño fácil y 
atractivo para el estudiante. Cada libro contiene objetivos de apren- 
dizaje muy explícitos, y cuestionarios de autoevaluación, así como 
respuestas y comentarios a la autoevaluación, permitiendo así al es- 
tudiante, darse cuenta de su propio progreso. Todos estos libros pue- 
den utilizarse como textos —ya sea en forma individual, o uniendo 
varias unidades— o como material de referencia. 


PELICULAS: Las películas han sido preparadas como material 
suplementario y de refuerzo para permitir al profesor mayor libertad 
de trabajo con cada estudiante. Las películas están en 16 mm, y en 
su mayor parte son en blanco y negro. A la fecha de publicación de 
estos libros sólo están a la venta en idioma inglés. 


CINTAS MAGNETOFONICAS: Al igual que las películas, estas 
cintas fueron preparadas como material suplementario y para ser uti- 
lizadas por el estudiante en su propia casa. A la fecha de publicación 
de estos libros sólo están a la venta en idioma inglés. 


NOTA A LOS MAESTROS: se han mantenido todas las referen- 
cias a otros elementos del programa en los materiales impresos, por 
considerarse que no obstaculizan la utilización de los libros como 
textos independientes. Los profesores que no tengan acceso a los ma- 
teriales audio-visuales fácilmente podrán sustituirlos con sus pro- 
pios experimentos o explicaciones. 


Esta obra se imprimió en los Talleres Gráficos de Carvajal « Cía. 
de Cali, Colombia, en el mes de abril de 1974. Se imprimieron 5000 
ejemplares en papel offset de 90 g/m?. 


4.1 


4.1.1 
4.1.2 
4.1.3 
4.1.4 


4.2 


4.2.0 
4.2.1 
4.2.2 
4.2.3 


4.3 


4.3.0 
4.3.1 
4.3.2 
4.3.3 
4.3.4 


4.4 


4.4.0 
4.4.1 
4.4.2 
4.4.3 


Tabla de materias 


Objetivos 

Diagrama estructural 
Glosario 

Notación 

Bibliografía 
Introducción 


Interpolación lineal y extrapolación 


La necesidad de una hipótesis 
Interpolación lineal 
Extrapolación 

Resumen 


Tablas de diferencia y operadores 


Introducción 

El operador diferencia 
Diferencias de polinomios 
Resumen 


Interpolación no lineal 


Introducción 

Ajuste de una función lineal 
Interpolación polinómica 
Fórmulas a partir de tablas 
Resumen 


Errores 


Introducción 
Equivocaciones 
Errores inherentes 
Resumen 


iii 


MB4.0 


OBSERVACIONES AL ALUMNO 

Al comienzo de cada texto encontrará lo siguiente. 

1. UNA PRESENTACION DE NUESTROS OBJETIVOS en la unidad, la cual puede usar para comprobar si ha completado el 
curso satisfactoriamente. 

2. UN DIAGRAMA ESTRUCTURAL que le muestra cómo están relacionadas las diferentes secciones de la unidad. El diagrama 
consta de: 
a) cuadros rojos que muestran la secuencia principal. 
b) cuadros negros punteados (a la izquierda del rojo) que le muestran lo que esperamos ya sepa antes de fijarse un estudio 

detallado. 

Cc) cuadros negros (a la derecha del rojo) que le ayudan a escoger temas que puede posponer u omitir si tiene poco tiempo. 

3. UN GLOSARIO de las definiciones, notaciones, etc., usadas en la unidad como referencia inmediata. 

4. UNA BIBLIOGRAFIA de los libros que puede consultar cuando haya acabado el trabajo de la unidad. 

5. CONCLUSIONES Y RESUMEN del trabajo realizado al final de cada texto. 


NUMERACION E INDICADORES 

Dentro de cualquier sección nos referimos a los temas por medio de un número. Al referirnos al mismo tema en otra sección, citamos 
cuatro números: número de la unidad, número del capítulo o parte, número de sección y número del tema. Por ejemplo, “Ecuación 3” 
si aparece en la misma sección, y “Ecuación. 4.3.1.3”, si está en otra sección. Los tres primeros números se citan en la parte superior 
de la página. Creemos que este es un sistema claro y fácil para consultar cualquier tema. Los indicadores serán una guía valiosa 
y le ahorrarán tiempo. 

En la margen derecha del texto hemos tratado de indicar la clase del material. Los términos usados son tan claros que no necesitan 
explicación. Por ejemplo, “Discusión” significa que una vez que haya entendido el tema, no necesita de nuevo recurrir a esta parte 
del texto. Estos INDICADORES tienen por objeto ayudarlo a localizar rápidamente las partes importantes cuando esté repasando 
o tenga que omitir parte del texto por falta de tiempo. 

También hemos tratado de indicar la importancia de los temas por el sistema de estrellas. 

*** Esto es muy importante. Usted debe entenderlo claramente. En algunos casos lo usará con tanta frecuencia, que se familiari- 
zará con ello sin necesidad de memorizarlo deliberadamente. El punto principal de este curso es no aprenderse los datos de me- 
moria; nosotros creemos que es más importante que entienda el trabajo y sea capaz de aplicarlo en la solución de los problemas. 
** Este es de menor importancia, pero tiene material al cual nos referiremos más adelante. 

» Esto generalmente no lo guía a trabajos adicionales en este curso y lo puede omitir si tiene poco tiempo. 

EW Indica el final de cada ejemplo, ejercicio o solución. 


Objetivos 


El objetivo general de esta unidad es describir algunos métodos prác- 
ticos de manejar datos en forma tabular. 


Después de estudiar esta unidad debe ser capaz de: 


(i) formar tablas de diferencias y usarlas para interpolación y ex- 
trapolación; 

(ii) explicar lo que significa hacer la diferencia de una función y po- 
der hacerla con algunas funciones elementales; 

(iii) usar métodos de diferencia para determinar y corregir equivoca- 
ciones en una tabla de valores; discutir los efectos de errores in- 
herentes en tablas de diferencia; juzgar cuántas diferencias pue- 
den emplearse útilmente en el método de Gregory-Newton; 

(iv) explicar el principio implicado en la interpolación polinómica; 

(v) usar la fórmula de Gregory-Newton para interpolación numérica 
(normalmente se espera que trabaje más allá de la aproximación 
cuadrática); 

(vi) usar métodos de diferencia para sumar series finitas de la forma 


> (a + br + cr?) 


r=1 


donde r toma valores enteros 
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Una función g se ajusta a una función f en 
los puntos x1,...Xn Si glxr) = flxx), pa- 
ak = 1 2... 


Las constantes de un polinomio se llaman 
coeficientes. 


El término continuo, aplicado a una función 
en este texto, significa que el gráfico de la 
función es continuo en el sentido de que no 
hay rompimientos o saltos en la curva. (En 
la Unidad 7, Sucesiones y límites I, se da 
una definición más rigurosa de continuidad.) 


Un cuártico es un polinomio de grado 4. 


El desarrollo binomial es la identidad poli- 
nómica 

n n 
(x+h)” =x" + [e + [Ejea Ho + 


(Véase RB9.) 


Si x1,...Xxm son puntos tabulares con se- 
paración h, igualmente espaciados, de una 
función tabulada f y si m > n, entonces las 
diferencias de orden n son las imágenes de 
X1, X2,-..Xm-=n bajo la función Af don- 
de A; es el operador diferencia de orden n 
para la separación h. 


Diferencias finitas es un término general 
aplicado a las primeras diferencias, segun- 
das diferencias, etc., y también a las diferen- 
cias flxa) — fíx1), flxa) — flxz),..., donde 
X1, X2,... están en orden ascendente pero 
no necesariamente a espacios iguales. 


Si Xm,...xn son puntos tabulares de una 
función tabulada f, entonces la extrapolación 
es el estimativo de las imágenes bajo f de 
los puntos que están fuera del intervalo 
[Xm, xn], por medio de una función g cons- 
truida que AJUSTA f a los puntos Xm,... Xn. 


La extrapolación lineal es extrapolación por 
medio de una función lineal. 


La extrapolación polinómica es extrapolación 
por medio de una función polinómica. 


La fórmula de Gregory-Newton para interpo- 
lación polinómica de grado n a la función f es 


8(x) = f(x) + OA, f(x) 
+30(0 — DARÍ(X) ++ 


+00 — 1)(0—2)...(0—n + DALf(x) 


30 


31 


34 


33 
33 


24 


40, 42 


FORMULA DE 
NEWTON-GREGORY 
FUNCION 
CUADRATICA 


FUNCION CUBICA 


FUNCION LINEAL 


FUNCION 
PERIODICA 


FUNCION 
POLINOMICA 


FUNCION REAL 


FUNCION TABULADA 


INTERPOLACION 


INTERPOLACION 
LINEAL 


INTERPOLACION 
POLINOMICA 


INTERPOLACION 
POLINOMICA 
DE LAGRANGE 


donde xx, Xx+1---Xisn son puntos tabu- 
lares igualmente espaciados con separación 


x—X 
Ly d= == Msbiel 
Xks 1 Xz 


Véase fórmula de Gregory-Newton. 


Una función cuadrática es una función poli- 
nómica de grado 2. 


Una función cúbica es una función polinó- 
mica de grado 3. 


Una función lineal es una función polinómi- 
ca de grado 1; ¡i.e. una función cuyo gráfico 
es una recta. 


Una función f con dominio R se dice perió- 
dica con período h si fíx + h) = f(x) para 
todo x E R donde h E R es el mínimo nú- 
mero positivo para el cual vale esta ecuación. 


Una función polinómica de grado n es una 
función de la forma x-——f(x), donde f(x) es 
un polinomio de grado n en x y el dominio de 
fes R o un subconjunto de R. 


Una función real es una función cuyo domi- 
nio y codominio es R o subconjuntos de R. 


Una función tabulada es una función real 
de la cual se han tabulado ciertos valores. 


Si xXm,...Xxn son puntos tabulares de una 

función tabulada f, entonces la interpola- 

ción es el estimativo de las imágenes bajo f 

de ciertos puntos que están dentro del inter- 

valo [xm, xn] por medio de una función 

construida g que ajusta f a los puntos Xm, 
Maa 


La interpolación lineal es interpolación por 
medio de una función lineal. 


La interpolación polinómica es interpolación 
por medio de una función polinómica. 


Si xx, Xk+1 Xx+2 SON tres puntos tabula- 
res de una función tabulada f, entonces el si- 
guiente polinomio cuadrático se llama poli- 
nomio de interpolación de Lagrange y ajusta 
ÍA Xr Xx+19 Xh+2 


(a — Xx 0 — Xx 4 2) 
(a — 4) 00 — Xy. +2) 


A Pl. 1) 


Vea =D 2) 


9 fl) 


Qs — A+ 2 — 1 


q(x) = F(xy) 
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El operador diferencia para separación h, de- 
notado por A,, opera sobre el conjunto de 
las funciones reales F, como sigue: 


ANI) = Lo — f(x + h) — f(x] (SEF) 


El operador diferencia de orden n para la se- 
paración h, denotado por As, opera sobre 
el conjunto de todas las funciones reales y 
es la composición de n veces el operador di- 
ferencia A, consigo mismo, i.e. 


A, = A, o A, o--:9 A, 


El término 0Í fíxx + h) — fíxx) | en inter- 
polación lineal se llama parte proporcional 
para el valor 6. En algunas tablas se han 
tabulado también las partes proporcionales, 
usualmente para 9 = 0,1, 0,2,..., 0,9. 


Las partes proporcionales promedio son par- 
tes proporcionales promediadas sobre un 
grupo de puntos tabulares. 


Un polinomio de grado n en x es una expre- 
sión de la forma anx” + an-1x”71 + 

+ ao, donde x se considera como una va- 
riable, an,... ag como constantes y an * 0. 


Un polinomio cuadrático es un polinomio de 
grado 2. 


Un polinomio cúbico es un polinomio de gra- 
do 3. 


Si x1,... xn son puntos tabulares igual- 
mente espaciados, con separación h, para 
una función tabulada f, entonces las prime- 
ras diferencias son las imágenes de x,,... 
Xn-1, bajo la función Af, donde A, es 
el operador diferencia para la separación h. 


Si f es una función tabulada, cuyas imáge- 
nes han sido tabuladas en correspondencia 
a elementos x;,..., xn del dominio, enton- 
Ces X1,..., Xn se llaman puntos tabulares 
(o argumentos tabulares). 


Véase diferencias de orden n. 


Una tabla de diferencias es una tabla cuya 
primera columna contiene los valores tabu- 
lares de una función tabulada, y cuyas co- 
lumnas siguientes contienen las primeras 
diferencias, segundas diferencias, etc. Los 
puntos tabulares usualmente están igual- 
mente espaciados. 


Los valores tabulares son las imágenes de 
los puntos tabulares bajo una función tabu- 
lada. 
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Notación 


Los símbolos se presentan en el orden en que aparecen en el texto. 


R 

log x 
ab 
ac Á 

q 

senx 

cos Xx 
az=b 

(a, b] 


A, 


ANS) o A, 
E 


Hof 


A; 


El conjunto de los números reales. 

El logaritmo de x en la base 10 (véase RB 11). 

La imagen de a bajo la función f es b. 

El elemento a pertenece al conjunto A. 

El conjunto de los enteros positivos. 

La imagen de x bajo la función seno (véase RB2). 
La imagen de x bajo la función coseno (ver RB2). 
a es aproximadamente igual a b. 


El subconjunto de R que consta de todos los números 
comprendidos entre a y b, incluidos a y b. 


El operador diferencia para la separación h. 

La función que es imagen de la función f bajo A». 
El conjunto de las funciones reales. 

La composición de las funciones f, y f2 (véase Uni- 
dad 1). 

El operador diferencia de orden n para la separación h. 
Una operación binaria (véase Unidad 3). 

La función f:x-——0 (x E R). 

a es mayor que b. 

La imagen de x bajo la función tangente (véase RB2). 
La suma y r”, de las potencias m de los primeros n 


r=1 


enteros positivos. 
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4.0 INTRODUCCION 


El propósito principal de esta unidad es decir a usted cómo trabajar 
con una función para la cual se han calculado algunos elementos de 
su gráfica, o se han dado. El concepto de función se explicó en la Uni- 
dad 1, Funciones; aquí tratamos con funciones cuyo dominio y codo- 
minio son conjuntos de números reales (i.e. R o subconjuntos de R). 
En principio, dado cualquier x que pertenezca al dominio, es siempre 
posible hallar la imagen de x aplicando la regla para la cual se definió 
la función. En la práctica, sin embargo, puede no ser posible usar este 
método cada vez que se requiera una imagen bajo la función: la regla 
puede ser muy complicada para aplicarla convenientemente, o en al- 
gunos casos, puede no ser conocida. En tales casos puede ocurrir que 
las imágenes (o valores aproximados de ellas) de algunos números 
particulares x,, x2, x3,..., del dominio puedan hallarse y en este 
caso estas imágenes se pueden enumerar junto Con xXx, X2, X3,..., 
en la forma de una tabla. Tal tabla es un subconjunto de la gráfica 
de la función. Los ejemplos más familiares de tales tablas son las ta- 
blas de logaritmos y de funciones trigonométricas que aprendemos a 
usar en el colegio. 


x log x 

1,0 0,0000 
1,1 0,0414 
1,2 0,0792 
1,3 0,1139 
1,4 0,1461 
1,5 0,1761 


Con frecuencia nos enfrentamos al problema de hacer un estimativo 
de las imágenes de valores de x que no han sido tabulados: por ejem- 
plo, ¿podemos, si es necesario, obtener un valor exacto de log 1,05 ó de 
log 1,6 a partir de los valores dados en la tabla 1? Este es el tipo de 
problema con el cual tenemos que tratar. 


Para funciones tales como el logaritmo, seno o coseno, el uso de las ta- 
blas es una cuestión de conveniencia solamente: Con la ayuda de un 
computador podemos calcular el logaritmo o el seno o el coseno de un 
número sin usar tablas (cómo se hace esto se explicará en la Unidad 
14, Sucesiones y límites II). Muchas funciones, sin embargo, están 
formadas a partir de reglas mucho más complicadas que la que define 
el logaritmo: un ejemplo sería la función: 


Fuerza neta hacia arriba por 


Posición sobre un , ' 
'——> ¿unidad de área del ala en esa 
ala de aeroplano nO 
posición 


Calcular las imágenes bajo esta función para un ala de aeroplano da- 
da (bajo condiciones de vuelo dadas), aun en unas pocas docenas de 
posiciones, sobre el ala, sería un proyecto importante de computación, 
y entre más posiciones se tomaran, más costoso sería el proyecto. Es 
mucho más práctico, por consiguiente, tabular las imágenes bajo la 
función en posiciones no muy numerosas y usar los métodos de esta 
unidad para hacer estimativos de las imágenes en otras posiciones. 
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* * 


Tabla 1 


Así que la revolución del computador no ha hecho obsoletas las téc- 
nicas de trabajo con funciones tabuladas (i.e. funciones para las cua- 
les se ha especificado en una tabla un subconjunto de su gráfico); 
más bien se ha extendido grandemente el conjunto de funciones a las 
cuales es posible aplicarlas. 


Hay otras situaciones también en las cuales nos podemos encontrar 
con funciones babuladas. La función tabulada puede resultar de una 
sucesión de observaciones experimentales, como por ejemplo, en las 
tablas usadas por los ingenieros para diseñar máquinas de vapor. 
Aquí se presenta un extracto de una de tales tablas*: 


Presión Punto de ebullición del agua 
(1b £/pul?) (F) 
300 417,33 
350 431,72 
400 444,59 
450 456,28 


En este caso la función está definida por una relación física en lugar 
de matemática, como en el caso del logaritmo, pero el dominio y el co- 
dominio de la función son también subconjuntos de R, así que el mé- 
todo de usar la tabla es el mismo de antes. Un ingeniero que necesi- 
tara el punto de ebullición del agua a una presión de, digamos, 325 
lIbf£/pul?, podría hacer un estimativo de él a partir de la tabla Il usan- 
do los mismos métodos que se usan para hacer un estimativo de log 
1,05 a partir de la tabla 1. Este problema de hacer un estimativo de la 
imagen de un número que está entre dos de los números tabulados en 
el dominio se llama el problema de interpolación. Uno de los objetivos 
de esta unidad es mostrar métodos exactos y convenientes para ha- 
cer interpolación. 


Un tercer tipo de función que puede tabularse es el que usted habrá 
encontrado si ha realizado alguna vez una prueba de inteligencia. Un 
ejemplo sencillo de una pregunta en tal prueba es el siguiente: 


¿Cuál es el número siguiente de la sucesión 
nl 


Esta sucesión se puede considerar como una lista de valores de imá- 
genes de una función que podemos llamar f, donde 


f:n ——(n-ésimo número de la sucesión), (n E 29), 


donde Z* es el conjunto de los enteros positivos. Por consiguiente 
corresponde a la tabla 


n  n-ésimo número de la sucesión 


UL tn 
sy UU Muy — 


o 


*O. W. Eshbach, Handbook of Engineering Fundamentals, (John Wiley 1966). 
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Tabla Il 


Definición 1 
a y 


Sucesión (1) 


El problema planteado en esta prueba de inteligencia es descubrir el 
valor de f(5), o sea la siguiente entrada de la tabla. Este es un proble- 
ma de naturaleza similar al de interpolación, excepto que el elemento 
del dominio cuya imagen estamos tratando de estimar no está entre 
dos de los valores tabulados sino más allá del último que se ha tabu- 
lado. El problema de hacer un estimativo de las imágenes bajo la fun- 
ción, para elementos del dominio que están fuera del conjunto tabula- 
do se llama problema de extrapolación. Otro ejemplo de extrapolación 
es el problema de obtener log 1,6 a partir de la tabla 1. Un segundo ob- 
jetivo de esta unidad es enseñar al lector cómo extrapolar. 


Observando la pregunta de la prueba de inteligencia que se muestra 
en la sucesión (1), usted probablemente se da cuenta de que el primer 
número que falta en la sucesión es 9. ¿Cómo llegó a esta conclusión? 
Puede haber razonado que el segundo número de la sucesión es mayor 
en 2 que el primero, que el tercero es mayor en 2 que el segundo, que 
el cuarto es mayor en 2 que el tercero, y por consiguiente que el quinto 
es probablemente mayor en 2 que el cuarto. Si razonó de esta manera, 
ya ha captado la idea básica del método que da el título a esta uni- 
dad, a saber, el de las diferencias finitas. En este método concentra- 
mos la atención sobre las diferencias entre los números sucesivos de 
una sucesión de imágenes bajo una función. Estas diferencias se lla- 
man diferencias finitas para distinguirlas de otras diferencias con- 
sideradas en una unidad posterior. En esta unidad posterior, Unidad 
12, Diferenciación I, tomamos la diferencia de imágenes de números 
del dominio a los cuales se permite estar arbitrariamente cercanos, 
mientras que en el tema de las diferencias finitas, la diferencia entre 
ellos se mantiene fija. Por razón del ligamen con la diferenciación, su 
estudio de diferencias finitas en esta unidad le proporcionará una in- 
troducción útil para algunas de las ideas de Diferenciación I y las 
otras unidades de cálculo, evitando por ahora el aparato extraconcep- 
tual (el concepto de límite) que se necesita para tratar con cantidades 
que pueden ser arbitrariamente cercanas. 


Como los computadores digitales están bien equipados para las dife- 
rencias finitas, hay procesos que se realizaban usando el cálculo y que 
ahora, en algunos casos, se calculan con la ayuda de las diferencias 
finitas. 


4.1 INTERPOLACION LINEAL Y 
EXTRAPOLACION 


4.1.1 La necesidad de una hipótesis 


Todo método de interpolación y extrapolación depende de una hipóte- 
sis acerca de la naturaleza de la función que se ha tabulado. Por ejem- 
plo, en la prueba de inteligencia que hemos mencionado: 


“¿Cuál es el siguiente número de la sucesión 
E A 


la respuesta obvia es 9, ya que los números parecen haberse formado 
de acuerdo con la regla “para hallar el número siguiente, agréguese 
2”: en otras palabras forman una progresión aritmética (sucesión) 


a,a+d,a+2d,a + 3d,... 


MB4.0/4.1.1 


Definición 2 
* * 


Definición 3 


* * 


4.1 


4.1.1 


Discusión 
* * 


Sucesión (1) 


Sucesión (2) 


cona = ly d = 2. ¿Pero es realmente tan simple como eso? ¿Por 
qué la respuesta no podría ser 11? Los números dados en la sucesión 
(1) podrían ser los primeros de una sucesión más larga que continúa: 


1,337, 11, 13,15, 121,22, 25,37, 3... 


Mirando solamente los cuatro números 1, 3, 5, 7 no tenemos manera de 
decir si provienen de la sucesión aritmética (2) o de la sucesión más 
complicada (3). Se cree, tanto por parte de los gestores de las pruebas 
de inteligencia como por aquellos que tienen la mala suerte de tener 
que responderlos, que la escogencia “razonable” es la sucesión arit- 
mética, quizás sobre la base de que la regla que la determina se puede 
hallar a partir de los números dados 1, 3, 5, 7. Pero, ¿podemos estar 
seguros de que existe una sola escogencia “razonable”? 


Este problema tipifica una situación que hallaremos varias veces en 
esta unidad. Para hacer una interpolación o (como aquí) una extrapo- 
lación, es necesario primero hallar una función simple que tenga los 
mismos valores en las imágenes que la función tabulada para los ele- 
mentos tabulados en el dominio. En el ejemplo de la prueba de inteli- 
gencia hay una buena posibilidad de adivinar la única función que 
tenía en mente quien aplicó la prueba, pero en general, como el objeto 
de la interpolación es con frecuencia aproximar una función complica- 
da por medio de una más simple, no habrá una manera única de esco- 
ger la función simple. 


El siguiente ejemplo ayuda a clarificar la última frase. 
Ejemplo 1 
Una parte de una tabla para la función seno 

x——> sen x (xeR) 


se muestra a continuación. 


pd sen Xx 
—0,4  —0,39 
03 03 
=0% <Q 
=01. 010 

0 0 

0,1 0,10 
0,2 0,20 
0,3 0,30 
0,4 0,39 


(Los valores de x están dados en radianes y las imágenes correspon- 
dientes están dadas con dos cifras significativas.) 


Ahora, parecería razonablemente claro que para cualquier elemento 
del dominio restringido [—0,4, 0,4], la función x-——x tiene aproxi- 
madamente la misma imagen que x-——sen x. Por consiguiente, si 
estuviéramos interesados solamente en este dominio restringido para 
la función seno, podríamos con confianza aproximar la función com- 
plicada 


xH—>sen x (x E [—0,4, 0,4]) 
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Sucesión (3) 


Ejemplo 1 
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por medio de la función simple 
sd (x ES [ —0,4, 0,4)) 


Sin embargo, como dijimos antes, esta aproximación no es única: por 
ejemplo se halla que tanto 


x—0os [5 —x] (x E [—0,4, 0,4]) 


como 
x? 
a (x E [—0,4, 0,4]) 
son también aproximaciones a la aplicación seno con este dominio res- 
tringido (la primera resulta ser exacta por cuanto 


cos 
Z 


La no unicidad no es algo que podamos eliminar, como se puede apre- 
ciar fácilmente considerando la posibilidad de asignar significados 
precisos a “complicado”, “simple” y “aproximada”. ES] 


3-> = sen x, (x E R)). 


Como una consecuencia de la no unicidad, hay diferentes métodos de 
interpolación, que bien pueden dar diferentes resultados. No hay un 
criterio absoluto para escoger entre estos métodos, excepto que, cuan- 
do el resultado se va a aplicar, digamos, a ingeniería de construcción, 
el criterio es: ¿funciona? Todo lo que puede hacerse es usar cualquier 
información que tengamos con respecto a la fuente de los números 
tabulados y escoger la nueva función de la manera más simple que pa- 
rezca conducir a una buena representación aproximada de la función 
original. Para hacer una buena escogencia es importante poder hacer 
un estimativo del error en un método de interpolación y extrapolación. 
Más adelante en esta unidad discutiremos brevemente cómo puede 
hacerse esto. 


Hemos mencionado que es importante saber de dónde provienen los 
números tabulados. El uso más importante de esta información es sa- 
ber si la función tabulada es suficientemente regular para justificar 
la aplicación de métodos de interpolación o de extrapolación. En la in- 
troducción mencionamos varios ejemplos de tales funciones “'regula- 
res”: el logaritmo, dos funciones que aparecen en ingeniería y una su- 
cesión de números de una prueba de inteligencia. Hay, sin embargo, 
funciones cuyo comportamiento es mucho menos regular y para las 
cuales los métodos descritos en esta unidad no son adecuados normal - 
mente. Estas son las funciones en las cuales toma parte un elemento 
de mucho azar o impredecibilidad. Un ejemplo extremo es la sucesión 
de números que se obtiene lanzando un dado o haciendo girar una ru- 
leta un gran número de veces. Aquí (afortunadamente para los pro- 
pietarios de casinos) aun las más sofisticadas técnicas de extrapola- 
ción son inútiles para predecir el número siguiente de la sucesión. 


Quizás la palabra “predecir” como se usó arriba necesita explicación. 
Si lanzamos un dado no cargado, podemos tener la certeza de obtener 
uno de los resultados 1, 2, 3, 4, 5,66 y en un gran número de lanza- 
mientos, cualquier conjetura de repetición de uno de esos seis núme- 
ros sería correcta aproximadamente una vez entre seis. Pero esto no 
es lo mismo que decir que si los primeros siete lanzamientos dan la 
sucesión 


6,3,6,2,1,4,4 


entonces el siguiente dará un 5. No hay ley que gobierne esta suce- 
sión: cada uno de los seis números es igualmente posible de aparecer 
en el siguiente lanzamiento. Este y otros tipos similares de sucesión 
se pueden clasificar informalmente como “sucesiones estadísticas”. 
Los métodos de diferencias finitas para interpolación y extrapolación 
no pueden aplicarse a este tipo de sucesiones. Son inapropiados por- 
que para usar estos métodos tenemos que hacer alguna hipótesis con 
respecto a la regla de formación de la sucesión. Por la misma natura- 
leza del azar implicado en el lanzamiento de un dado, no podemos ha- 
cer este tipo de hipótesis para las sucesiones estadísticas. Volvere- 
mos a estudiar este tipo de sucesión en la Unidad 16, Probabilidad y 
estadística l. 


Las sucesiones estadísticas son de un tipo menos extremo en los da- 
tos económicos, tales como cifras comerciales por trimestres, o precios 
en la Bolsa. Por ejemplo, a continuación se da una tabla del precio de 
acciones en Consolidated Gold Fields durante la semana del 9 al 16 
de agosto de 1968: 


Día Precio Cambio 
Viernes 9 de agosto  66/6 
Lunes 12 - 68/- +1/6 
Martes 13 69/- +1/- 
Miércoles 14 70/- +1/- 
Jueves 15 72/- +2/- 
Viernes 16 74/- +2/- 


El arte de hacer dinero en la Bolsa es el arte de predecir el comporta- 
miento de una sucesión como la de la tabla; pero aunque los números 
de la tabla no provienen del azar como los de lanzar un dado, es poco 
posible predecir la sucesión por métodos de diferencia finita. Los nú- 
meros de la columna de “cambio” (lo que llamamos “diferencias” en 
esta unidad) son todos por lo menos 1/—, lo que podría llevarlo a uno 
a esperar que los números continúen creciendo firmemente, quizás en 
unos 5 ó 10 chelines en la semana siguiente. En efecto, al final de la 
semana siguiente el precio fue el mismo que en la 16* y al final de la 
semana siguiente regresó a 69/—. Lamentablemente, no podemos 
prometerle que el estudio de esta unidad le ayudará a hacer fortuna 
en la Bolsa. La razón es que los precios de las acciones están deter- 
minados por muchos otros factores además de los que se muestran en 
la tabla. Desde luego hay “tendencias”, pero éstas sólo pueden iden- 
tificarse confiablemente por otros métodos; volveremos a este tema 
en las unidades sobre probabilidad y estadística. 


La conclusión principal que puede sacarse de la discusión de esta sec- 
ción es la de que los métodos de interpolación y extrapolación traba- 
jan mejor cuando la tabla de números que se consideren provenga de 
una función matemática bien definida o empírica de la cual pueda 
esperarse que conduzca a un patrón aceptablemente regular de varia- 
ción en los valores imagen. Comenzamos por observar algunos casos 
en los cuales este patrón es suficientemente simple para ser aproxi- 
mado por medio de una función lineal, esto es, una cuyo gráfico es una 
recta. Se ha considerado ya un ejemplo que supusimos ser de este ti- 
po: la extrapolación de la sucesión 1, 3, 5, 7,... para la cual la función 
aparente 


k—— (k-ésimo número de la sucesión) 
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Definición 1 
* * 


tiene el siguiente gráfico: 


k-ésimo número 
de la sucesión 


9 *— Extrapolado 


puntos dados 


12345 k 


En la siguiente sección consideraremos el método de interpolación co- 
rrespondiente. 


Ejercicio 1 


Escribir un párrafo de discusión sobre la aplicabilidad de los métodos 
de diferencias finitas a los siguientes datos tabulados. 


Elemento químico Número atómico Peso atómico 
Manganeso 25 54,93 
Hierro 26 55,85 
Cobalto 27 58,94 
Níquel 28 58,69 
Cobre 29 63,54 
Zinc 30 65,38 
Tono de nota musical Frecuencia de vibración* 

Do mediano 256 

Re 288 

Mi 320 

Fa 3413 

Sol 384 

La 4263 

Si 4903 

Do 512 


Miles de hombres en armas: Ejército, marina y fuerza aéreat 


A A A A A 
Año Inglaterra Francia Rusia E.E.U.U. Alemania Italia 


e. /K/K/<Ká 


1914 397 834 1254 165 864 345 
1924 349 835 678 251 114 359 
1932 316 690 562 237 114 324 
1939 460 864 2269 334 1182 1077 
1949 770 589 4000 1616 — És 
1955 803 950 4500 2935 — E 


HA E A 


«Sir James Jeans, Science and Music, (Cambridge University Press 1961). 
+P. J. Noel-Baker, The Arms Race, (Oceana Publications 1956). 
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Ejercicio 1 
(10 minutos) 


Solución 1 


Para cada una de las tres tablas su respuesta debe considerar los si- 
guientes aspectos: 


(i) Si están determinados los valores tabulados (elementos del codo- 
minio de alguna función no especificada) por los puntos tabulares 
(elementos del dominio de esta función) de acuerdo con alguna 
ley matemática, científica o de otro tipo natural o si están afec- 
tados significativamente por otros factores (incluido el azar). 

(ii) Si la respuesta a (i) es que están determinados por alguna ley, 
debe determinarse si esta ley es suficientemente simple y regular 
para aproximarla por medio de una expresión matemática simple 
obtenida de los números de la tabla. 


En el ejemplo de química, hay sin duda una ley científica, pero no es 
lo suficientemente regular para justificar el uso de los métodos de di- 
ferencias finitas. (Obsérvese cómo los pesos atómicos en su mayor 
parte crecen con el número atómico, pero esto falla entre el cobalto y 
el níquel). En el ejemplo de música las frecuencias de vibración están 
determinadas por las notas de acuerdo con leyes acústicas, y la ley es 
lo suficientemente regular como para justificar el uso de las diferen- 
cias finitas; en efecto, desde la invención de la escala igualmente 
temperada (que es una buena aproximación de la escala diatónica 
dada en el ejercicio) el método usual de afinar un piano se ha basado 
en una sencilla fórmula matemática (cada semitono corresponde a un 


factor de J2 6 sea 1,059 en frecuencia). En el ejemplo de la milicia, 
la magnitud de la fuerza armada de un país sin duda influenciada 
por factores distintos del simple paso del tiempo y de las magnitudes 
de las fuerzas armadas de los otros países; así que la tabla, (como la 
de los precios de las acciones dada en el texto) no proporciona infor- 
mación suficiente para justificar el uso de métodos de diferencias fi- 
nitas; la pregunta (ii) no se aplica en este caso. ls] 


4.1.2 Interpolación lineal 


El método más simple de interpolación es el que se usa comúnmente 
en las tablas de funciones elementales, tales como las tablas de lo- 
garitmos. En la introducción de esta unidad, por ejemplo, menciona- 
mos el problema de hallar el log 1,05 a partir de la tabla siguiente: 


Xx log x 
1,0 0,0000 
1,1 0,0414 
1,2 0,0792 
1,3 0,1139 
1,4 0,1461 
1,5 0,1761 
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Solución 1 


4.1.2 


Discusión 
* * 


Tabla 1 


Como 1,05 está a mitad de camino entre 1,00 y 1,10, el estimativo más 
simple es suponer que log 1,05 está a mitad de camino entre log 1,00 y 
log 1,10, o sea 


0,0000 + 0,0414 


log 1,05 = 
og 1,0 3 


= 0,0207 


Esto está de acuerdo en 3 cifras decimales con el valor 0,0212 dado en 
las tablas de logaritmos con cuatro decimales. (Este último valor es a 
su vez correcto hasta cuatro cifras decimales; el valor exacto de log 
1,05 es un decimal que no termina). De manera similar, para hallar 
log 1,13, puesto que 1,13 está a ¿5 del camino entre 1,10 y 1,20, el esti- 
* mativo más simple es log 1,10 más $; de la distancia entre log 1,10 y 
log 1,20, esto es 


log 1,13 = log 1,10 + ¿(log 1,20 — log 1,10) 


0,0414 + +$5(0,0792 — 0,0414) 
= 0,0527 


Esto está de acuerdo hasta tres cifras decimales con el valor 0,0531 
dado en las tablas de cuatro cifras. 


El método usado para deducir las ecuaciones (1) y (2) se llama a ve- 
ces método de partes proporcionales, ya que el número de la derecha 
de la ecuación (2) divide el intervalo [log 1,10, log 1,20] en la mis- 
ma proporción que el número 1,13 divide el intervalo (1,10, 1,20]. 


x log 
1,00 0,0000 
1,05 A 
1,10 0,0414 
1,13 QQ 
1,20 0,792 


Para tratar con el caso general, denotemos la función tabulada por f y 
los números del dominio tabulado por x;, X2,... Xn, Ordenados en 
forma creciente. Estos se llaman puntos tabulares. Las imágenes de 
estos puntos son entonces fíx1), flx2),..., f(xn); se llaman valores 
tabulares de la función. Por brevedad denotaremos los valores tabula- 
Tes Por Y1, Y2,---> Yn! 


x J(x) 


X1 yy =f(x1) 
Xx Ya = f(x,) 


Xn Ya = f(x) 
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Ecuación (1) 


Ecuación (2) 


Tema principal 
h $ £ 


Definición 1 
h * * 
Definición 2 
+ * $ 
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Para hacer un estimativo de f(x) cuando x está entre dos puntos ta- 
bulares x, y x2, digamos, hallamos la proporción en la cual el núme- 
ro x divide al intervalo [x,, x2] y hallamos el número y que divi- 
da al intervalo [y,, y2] en la misma proporción. Podemos repre- 
sentar el procedimiento gráficamente: 


B_ Gráfico de f 


La proporción en la cual el número x divide 'al intervalo [x;, x>2] 
es 


x—=X, _AD 
Xa -X, AC 


Usando el hecho de que los triángulos ADE y ACB son semejantes, 
tenemos 


Como CB es la longitud del intervalo [y,, y21, E es el punto cuya 

coordenada y divide al intervalo [y;,, y2] en la proporción reque- 

rida. Así vemos que nuestro método de hacer un estimativo de f(x) en- 

tre A y Bes aproximar su gráfica por medio de una recta entre esos 

dos puntos. Por esta razón este método se llama de interpolación li- Definición 3 
neal. o 


Algunas tablas impresas dan alguna ayuda sobre el método de inter- 
polación lineal colocando los valores de yx+, — yx a la derecha de 
los valores de y» y yx+1 y sobre una línea a mitad de camino entre 
ellos, como por ejemplo en el siguiente extracto de una tabla de recí- 
procos: 


A _ A _ _—— 


Xx 1/x 
1,60 0,6250 > 
1,61. 0,6211 7 se 
1,62 0,6173 7 38 
1,63. 0,6135  — 


áK A C—6 


Los números —39, —38, —38, que se muestran como primeras diferen- Definición 4 
cias, deberían ser en realidad —0,0039, — 0,0038, —0,0038, pero usual- *.o ros 
mente aparecen como se ha expresado (o aún sin el signo menos) 

para ahorrar espacio. 
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Ejercicio 1 
Usar interpolación lineal para obtener valores aproximados para 


(i) log 1,25, usando la tabla I de la Sección 4.1.2; 
(ii) el punto de ebullición del agua a una presión de 365 lbf/pul? 
a partir de la siguiente tabla 


Presión Punto de ebullición del agua 


(Ib£/pul?) (*F) 
300 417,33 
350 431,72 
400 444,59 
450 456,28 


El método de interpolación lineal puede aplicarse a cualquier tabla, 
pero sólo en casos especiales se puede confiar en los valores así obte- 
nidos. El criterio es que el gráfico de la función debe ser aproximada- 
mente una recta entre los puntos tabulares. 


Ejercicio 2 


Tabular los valores de sen x para x = 0, rr, 21, 3r radianes y 
T % y % , 

calcular un valor de sen — por interpolación lineal. ¿Qué se puede con- 

cluir? 10 El 


Ejercicio 3 


(i) Tabular los valores de sen x para x = 0, a y r, y calcular 
un valor para sen ¿6 Por interpolación lineal. 
1U 
(ii) Calcular un valor para sen . por 

. . De . . 10 Xx sen x 

interpolación lineal a partir de la 

tabla de la derecha. 0 0 
T 
— 1 
16 0,1951 
2x1 
— 0,3827 
16 
3 
e 0.5 
16 0999 
58 0,7071 
16 


(iii) Con cuatro cifras decimales, el valor de sen 0 es 0,3090. ¿Qué 


concluye usted comparando esto con los resultados del ejercicio 
anterior y (i) y (ii) de este ejercicio? 5 
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Ejercicio 1 
(3 minutos) 


Discusión 
Rh * 


Ejercicio 2 
(3 minutos) 


Ejercicio 3 
(3 minutos) 


Solución 1 


(i) 0,0966 
(ii) 435,58 *F. 


Solución 2 


X 0 T 2x1 3 
sen x 0 0 0 0 
La interpolación lineal da sen Ñ = (). En 


lineal da sen x = 0 para todos los valores 
3”. La conclusión es que la interpolación 


realidad, la interpolación 


x comprendidos entre 0 y 
lineal puede dar resulta- 


dos sumamente inexactos si el intervalo de tabulación es demasiado 
grande. Se puede concluir también que debe tenerse mucho cuidado 
en los casos en que se sabe poco acerca de la función tabulada. MW 


Solución 3 


(i) 


Xx sen Xx 

0 0,0000 
A 0,7071 
S 1,0000 
= 0,7071 
- 0,0000 


El valor interpolado para x = El es 0,0000 + (0,7071 — 0,0000)0, con 


10 
T 
mm % , 
As ==, 
0 10 
4 


= 0,0000 + 0,7071 x a = 0,2828 
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(ii) El valor interpolado para x = q es 


T T T 
sen E + [sen[2=2 sen) o 
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Solución 1 


Solución 2 


Solución 3 


con 
E: 

10 16_6 

72 .» 10 
16 16 


= 0,1951 + [0,1876 x $) = 0,3077 


(iii) Los valores interpolados se hacen más cercanos al valor de cua- 
tro cifras decimales cuando reducimos el intervalo entre los pun- 
tos tabulares. E 


El ejercicio 2 ilustra un caso en que la interpolación lineal es comple- 
tamente inapropiada: la función seno no puede ser representada sa- 
tisfactoriamente por una recta que una puntos tabulares tan lejanos 
como Ú y r, como es obvio en la figura. 


sen x 


Gráfico de x +» sen x 


t 
interpolación 
lineal 


Tr 
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Por otra parte, el ejercicio 3 muestra dos casos en que los puntos tabu- 
lares son más cercanos y donde la interpolación lineal da resultados 


aceptables. En (i) la separación entre los puntos tabulares es 3 y la 
interpolación tiene un error de 0,3090 — 0,2828 = 0,0262; en (ii) la 
anchura del intervalo es sólo de = y la interpolación tiene un error 
de solamente 0,3090 — 0,3077 = 0,0013. Así que, en este ejemplo, la 


exactitud de la interpolación lineal es generalmente mejor, mientras 
más cercanos están los puntos tabulares. 


sen x 


Gráfico de x ==» sen x * 
interpolación con intervalo =- 


interpolación con intervalo > 
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Discusión 
* 


Al diseñar tablas matemáticas tales como las de logaritmos, usual- 
mente se escoge la separación tabular suficientemente pequeña para 
que el error de la interpolación lineal no sea mayor que el de redondeo 
en las cifras tabuladas. Esto implica que una mayor exactitud en las 
tablas exige una separación tabular más pequeña y por consiguiente 
más entradas tabulares. Por ejemplo en Chambers's Four-Figure Ta- 
bles (Tablas de Chambers con cuatro cifras decimales) los logaritmos 
ocupan 4 páginas, pero en Chambers's Six-Figure Tables (Tablas de 
Chambers con seis cifras decimales) ocupan 26 páginas. En las tablas 
con cuatro cifras se acostumbra a hacer la interpolación más conve- 
niente, sacrificando en algo la exactitud, remplazando las partes pro- 
porcionales por partes proporcionales promedio, que son simplemente 
las partes proporcionales promediadas sobre un grupo de intervalos 
tabulares (en lugar de tomar los puntos tabulares más cercanos). Por 
ejemplo, en la tabla de recíprocos de esta página, la parte proporcional 
media bajo el número 2 de la tabla se obtiene así: 


(5882 — 6250) Ñ 2 _ m6. AÑ 
10 10. 1007 


(El signo menos no se muestra en la tabla por brevedad; se espera que 
1 ' 
usted recuerde sustraer, ya que .s decrece cuando x crece.) 


Así que mientras está leyendo el recíproco de 16,12 ó de 16,72 usa la 
misma parte proporcional media, i.e. —7. Pero si usara los puntos ta- 
bulados más próximos debería usar interpolación lineal sobre los va- 
lores tabulados en 16,1 y 16,2 para 16,12, para obtener la verdadera 


parte proporcional, 0,2(6173 — 6211) = —8, en vez de —7 y sobre los 
valores tabulados en 16,7 y 16,8 para 16,72 para obtener la verdadera 
parte proporcional 0,2(5952 — 5988) = —7. 

Ejercicio 4 


Usando la tabla de recíprocos, leer los recíprocos de 1,66 y 1,67 (en 
tales tablas se espera que usted halle la posición del punto decimal 
por sí mismo). Por interpolación lineal, obtener un valor para el recí- 
proco de 1,667. Obtener también un valor del recíproco usando las 
partes proporcionales promedio dadas y explicar cualquier discre- 
pancia. Bl 


Tabla de recíprocos 


16 | 6250 -6211 6173 6135 6098 6061 6024 5988 5952 5917 


Partes proporcionales promedio 


12 34 5 6.-8,.8 9 


4 7 11 15 18 22 26 30 33 
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Ejercicio 4 
(5 minutos) 


4.1.3 Extrapolación 


Lo más importante que debe saberse con respecto a la extrapolación es 
que puede ser considerablemente menos exacta que la interpolación. 
Para la extrapolación lineal esto puede verse en la figura siguiente, 
que muestra el gráfico de una función junto con la recta que la cruza 
en dos puntos tabulares xs y xx+1. Para valores de x en el interva- 
lo [xx, xx+1], esta recta corresponde a la fórmula de interpola- 
ción dada en la sección anterior y puede ser más bien exacta si el in- 
tervalo no es demasiado ancho. Sin embargo, si la recta se usa para 
extrapolación, que está fuera de este intervalo, la exactitud disminu- 
ye rápidamente. 


t() 


Error más 
pequeño 


Aproximación con 
una recta 


Y rc 


Variación de la Variación de la Xx+1 Variación de la 
extrapolación interpolación extrapolación 


| 
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Por esta razón, si hay una escogencia entre interpolación y extrapola- 
ción, es preferible siempre la interpolación. A veces, sin embargo, no 
podemos evitar la extrapolación. Esto puede ocurrir, por ejemplo, si 
estamos calculando entradas sucesivas en una tabla de alguna fun- 
ción nueva y necesitamos un estimativo más o menos aproximado de 
la siguiente entrada de la tabla para usarlo como primera aproxima- 
ción en un proceso iterativo para calcularla exactamente. Este es el 
problema que consideraremos aquí: dada una tabla incompleta de una 
función, hacer un estimativo de la siguiente entrada de la tabla. 


Hemos considerado ya un problema de este tipo: la pregunta de la 
prueba de inteligencia en la cual había que adivinar el siguiente nú- 
mero de la sucesión 1, 3, 5, 7,... 


El método por el cual resolvió este problema se puede sistematizar ta- 
bulando la sucesión y anotando las diferencias de los valores tabula- 
res sucesivos como en las tablas impresas que se dieron como ejemplo 
para la interpolación lineal: 


k k-ésimo número 
de la sucesión 


UI hh Ynun — 
Y Mw 
mM inn 


15 


MB4.1.3 
4.1.3 


Tema principal 
R *$* * 


(continúa en la página 16) 


MB 4.1.2/4.1.3 
Solución 4.1.2.4 Solución 4.1.2.4 
1/1,66 = 0,6024, 1/1,67 = 0,5988 
Interpolación lineal 
1/1,667 = 0,6024 + (0,5988 — 0,6024)0 
con 


_ 167 —- 1,66 7 
- 1,667- 1,66 10 


- 
de modo que 1/1,667 = 0,6024 — 0,0036 x 10 = 0,5999 


Partes proporcionales promedio 
1/1,667 = 0,6024 — 0,0026 = 0,5998 


(La parte proporcional promedio para 0 = ¿5 se muestra en la tabla 
como 26). 


Los resultados que se obtienen usando partes proporcionales prome- 
dio son ligeramente inexactos porque las partes proporcionales para 
15 varían entre 


(6250 — 6211) x %; = 27 


(5952 — 5917) x $ = 24 


a través de la tabla, de modo que usando el valor medio 26 podemos 
introducir un error hasta de dos dígitos en el último puesto decimal. Mi 


(viene de la página 15) 


La columna de doses está integrada por las diferencias de números 

sucesivos de la sucesión. La manera natural de continuar esta colum- 

na es insertar otro ,; esto implica que la entrada siguiente en la su- 

cesión debe ser un y. Este es un ejemplo de extrapolación lineal» UN Definición 1 
proceso que puede usarse siempre que las entradas de la columna de e e é 
diferencias son todas iguales. 


Ejercicio 1 Ejercicio 1 
: e e : ó (2 minutos) 
Usar interpolación lineal para hallar estimativos de los valores tabu- 


lares que faltan en la tabla siguiente. 


AS sen x 
29" 54' 

30? 0,5000 
30* 6' 0,5015 
30* 12 0,5030 
30? 18' 0,5045 
30* 24' 


Un problema de extrapolación menos trivial es el que da la tabla de 
distancias de parada impresas en la contratapa del “Código de carre- 
teras”: 


Rapidez Distancia de parada e 
(millas/h) (pies) 
20 40 
30 75 
40 120 
50 175 
60 240 


Supongamos que el examinador de conducción de vehículos le pide la 
distancia de parada para 70 millas/h o para 10 millas/h. Aquí hay 
una manera de hacer un estimativo de ellas. Formamos otra vez la co- 
lumna de diferencias 


40 
35 

75 
45 

120 
55 
175 6s 

240 


Esta vez las diferencias no son iguales, pero es más bien fácil adivi- 
nar la siguiente entrada de la columna de diferencias y completar así 
la solución. 


Ejercicio 2 
¿Cuál es la distancia de parada para una rapidez de 70 mi/h? El 
(La respuesta a este ejercicio está contenida en lo que sigue.) 


Las diferencias de la tabla II forman una sucesión aritmética con di- 
ferencia común 10; por consiguiente la manera natural de continuar 
la columna de diferencias es colocar 75 en la entrada siguiente. En- 
tonces la entrada siguiente en la columna de distancia de parada de- 
be ser 315, ya que 315 — 240 = 75. El proceso por el cual obtuvimos 
75 en la entrada siguiente de la columna de diferencias se puede sis- 
tematizar construyendo otra columna de números, las diferencias de 
las entradas sucesivas en la columna de diferencias. Estas se llaman 
segundas diferencias de la sucesión original de distancias de parada. 
Con las segundas diferencias incluidas la tabla queda: 


Distancia Primeras Segundas 
Rapidez de parada diferencias diferencias 


20 40 
30 75 ie 10 
40 120 po 10 
50 175 ió 10 
60 240 3 10 
70 315 
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Tabla Il 


Ejercicio 2 
(2 minutos) 


Tema principal 
h «* 


Tabla III 


(continúa en la página 18) 


Solución 1 


La extrapolación lineal da sen (30? 24') = 0,5060 y sen (29” 54') 
0,4985. 


(viene de la página 17) 


Como las segundas diferencias de los números dados son todas 10, es- 
peramos que la entrada siguiente en la segunda columna de diferen- 
cias sea también 10 y las otras cifras que aparecen en rojo se determi- 
nan entonces resolviendo las ecuaciones 


15-65=10 y  315-240=75 


El arreglo de números de las últimas tres columnas de la tabla III se 
llama una tabla de diferencias. Ml 


Ejercicio 3 
Use las cifras del “Código de carreteras” para hacer un estimativo de 
las distancias de parada para una rapidez de 80, 10 y 0 mi/h. SE 
Ejercicio 4 


He aquí una parte de la tabla de cuadrados. 
Formar la tabla de diferencias para ella y ha- x Xx 
cer un estimativo de 37? por extrapolación. 


31 961 
32 1024 
33 1089 
34 1156 
35 1225 
36 1296 
37 As 
E 


Ejercicio 3 
Las imágenes bajo una función f están tabuladas para valores igual- 
mente espaciados de la variable. Si las segundas diferencias son 


(i) todas cero, 
(ii) no todas cero, 


¿qué se puede deducir acerca del gráfico de la función? A 


Ejercicio 6 (Omiítalo si hizo el ejercicio 5 completamente bien). 


Tabule las imágenes bajo la función 

x——>x +senrx (xeR) 
para x = 0, 1, 2, 3 y calcule las segundas diferencias de esta tabla. 
Calcule las imágenes también para x = 3,3,3 y luego trace la gráfica 


de la función. Este es un ejemplo de una tabla con segundas dife- 
rencias cero a partir de una función que no es lineal. E 
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Solución 1 


Definición 2 
EY 
Ejercicio 3 
(2 minutos) 


Ejercicio 4 
(2 minutos) 


Ejercicio 5 
(5 minutos) 


Ejercicio 6 
(3 minutos) 


4.14 Resumen 


En esta parte hemos mostrado que las imágenes de una función tabu- 
lada correspondientes a valores comprendidos entre valores tabula- 
res conocidos, se pueden calcular convenientemente usando interpo- 
lación lineal, previsto que dentro del intervalo que se usa para inter- 
polación el gráfico de la función sea lo suficientemente próximo a una 
recta. Las imágenes correspondientes a valores situados fuera de los 
puntos tabulares también pueden hallarse de manera similar, pero la 
extrapolación puede ser peligrosa a menos que tengamos algún cono- 
cimiento de la función fuera del recorrido tabulado. 


4.2 TABLAS DE DIFERENCIA Y OPERADORES 


4.2.0 Introducción 


Vimos en la sección anterior que tanto para interpolación como para 
extrapolación es útil agregar una tabla de diferencias (o sea, una o 
más columnas de diferencias) a los valores tabulados de una sucesión 
o función. Para interpolación y extrapolación lineales la tabla de dife- 
rencias necesita incluir solamente las primeras diferencias de los 
valores tabulares de la sucesión o función, pero cuando los métodos 
lineales no son lo suficientemente exactos porque el gráfico de la fun- 
ción no es una línea recta (como en el ejemplo del “Código de carrete- 
ras” usado en la sección precedente; véase también ejercicio 4.1.3.5 
(i)), es posible que necesitemos incluir en la tabla columnas de se- 
gundas diferencias y quizás también de terceras diferencias (las dife- 
rencias de entradas sucesivas en la columna de segundas diferen- 
cias) o aun diferencias de orden superior. La utilidad de las tablas de 
diferencias no está, en manera alguna, confinada a la extrapolación. 
También las usamos para interpolación; los métodos resultantes pue- 
den tener mayor exactitud que el de interpolación lineal discutido an- 
tes. Ellas hacen esto posible porque las segundas diferencias dan una 
indicación de la desviación de la gráfica con respecto a la recta, por- 
que muestran cuán rápido van cambiando las primeras diferencias a 
medida que nos desplazamos sobre la tabla (véase la figura) y pueden 
usarse para corregir esta desviación. 


A 
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(continúa en la página 21) 


MB4.1.3 
Solución 4.1.3.3 Solución 4.1.3.3 
A partir de la tabla obtenemos los siguientes estimativos: 
a 80 mi/h la distancia de parada estimada es 400 pies. 
a 10 mi/h la distancia de parada estimada es 15 pies. 


a 0 mi/h la distancia de parada estimada es 0 pies. n 
Solución 4.1.3.4. Solución 4.1.3.4 
Primeras Segundas 
x y diferencias diferencias - 
31 961 
32 1024 ja 2 
33 1089 67 2 
34 1156 69 2 
35 1225 31 2 
36 1296 73 2 
37 
1369 s 
Solución 4.1.3.5 Solución 4.1.3.5 


(i) Los puntos de la gráfica correspondientes a puntos tabulares del 
dominio están sobre una recta, pero el gráfico mismo no es nece- 
sariamente una recta. 

(ii) El gráfico no es una recta. Es 


Solución 4.1.3.6 Solución 4.1.3.6 


Primeras Segundas 
x x+senrmx diferencias diferencias 


0 0 1 
1 1 1 0 
2 2 1 0 
3 3 
A 1,5 
3 0,5 
3 3,5 
t (0) 


Gráfica de x= x +sentrx (xe[0,3]) 
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"(viene de la página 19) 


Ys - Ys 


muy 
desiguales 


Y: - Y: 


«“— h — — h — Ñ *x 


La interpolación no lineal se considerará en la Sección 4.3 de esta 
unidad. 


Otra aplicación de la tabla de diferencias que algunas veces es útil, es 
detectar una equivocación en los valores tabulados de la función. Es- 
ta aplicación se considerará en la Sección 4.4 de esta unidad. 


Antes de que podamos tratar con cada una de estas aplicaciones, sin 
embargo, necesitamos observar la estructura de la tabla de diferen- 
cias misma con más detalle. 


4.2.1 El operador diferencia 4.2.1 


La estructura de la tabla de diferencias se determina por dos factores: Tema principal 
cómo están construidos los valores tabulares originales y cómo está *orsx 
construida la tabla de diferencias a partir de estos valores tabulares. 
Consideramos primero los valores tabulares originales. Estos son las 

imágenes de los puntos tabulares x,, x2,..., xn bajo la función da- 

da f. Hasta ahora los puntos tabulares han sido ordenados solamente 

por magnitud. Ahora, para simplificar la teoría, haremos explícita una 

condición que estaba en efecto satisfecha por todos los ejemplos que 

hemos considerado: que los puntos tabulares estén igualmente espa- 

ciados. La separación entre los puntos tabulares puede ser cualquier 

número positivo; se denota usualmente por h (otra manera de decir 

esto es que las primeras diferencias de la sucesión de puntos tabula- 

res sean todas iguales a h). Con esta simplificación, tanto los pun- 

tos como los valores tabulares se forman con aplicaciones definidas a 

los puntos tabulares y tenemos la estructura que se muestra a conti.- 

nuación, donde las flechas horizontales indican una aplicación de la 

función f y las verticales una aplicación de la función 


sumar h:x——>x+h (xeR) 
A 4 1(x2) 


E ES 
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El segundo factor que determina la estructura de la tabla de diferen- 
cias es el método de formar la tabla misma por repeticiones de una so- 
la operación binaria básica; esta operación consiste en sustraer dos 
entradas sucesivas en una columna y colocar el resultado a la derecha 
como se muestra en el siguiente diagrama. 


f(x.) 


a 1(x,) Xx 
A + a 
y 10) — f(x) 


1(5) 


Este diagrama da una descripción de la estructura de cualquier tabla 
de diferencias; la operación básica que usa es una operación binaria 
que actúa sobre dos números de la tabla. 


En el caso que estamos considerando aquí, con los puntos tabulares 
igualmente espaciados, hay otra manera de considerar las columnas 
de diferencias. Para hacer esto hacemos uso del hecho que x2 — x¡ =h, 
x3 — x2 = h, etc., de modo que el diagrama anterior se puede recons- 
truir así 


f(x.) 


a o 
¿7 A 
7 hal ES 1 + h) - f(x) 


f(x.) 


Esta representación muestra que la columna de las primeras diferen- 
cias de f se puede considerar como una tabla de la nueva función 


afín + — Noe der] 


con los mismos puntos tabulares x,, x2,... 


Aquí tenemos un operador del tipo discutido en la Sección 1.1.7 de 
la Unidad 1, Funciones. Dicho operador toma cualquier función f con 
dominio en R y la aplica a la nueva función definida antes. Como para 
h dado la nueva función está unívocamente determinada por f, el ope- 
rador es una función cuyo dominio es un conjunto de funciones: el 
conjunto de todas las funciones con dominio y codominio en R. Este 
operador se denotará por 3,, y se llama operador diferencia Para la 
separación h. La fórmula que define el operador diferencia es 


Ay: —[x > f(x + h) — f(x) (x,x + hedominio de f)| (feF) 
donde F es el conjunto de todas las funciones reales. (Definimos co- 
mo función real aquella cuyo dominio y codominio son conjuntos de 


números reales.) Si omitimos el dominio de la función imagen por 
brevedad, esta fórmula se puede escribir más brevemente así 


AN [x—> f(x + h)- fx) (SeF) 


159) 
159] 
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Definición 1 
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Notación 1 


Definición 2 
* * ; 


Notación 2 


Una forma alterna de formular la definición de A, es especificar la 
imagen de cualquier elemento de su dominio: 


AN = Lo f(x +) JN (SEF) 


donde también se omite el dominio de la función imagen por brevedad. 


Ejemplo 1 
Considérese 
f:x—1+2x+x? (xeR) 


Queremos hallar Ao, (f). Obsérvese que ahora hemos especificado 
h,i.e. h = 0,2. 


Por consiguiente 

Fx +h)= f(x + 0,2) = 1 + 2x + 0,2) + (x + 0,2)2) 
Se sigue que 

Ay f) = [lx —— (1 + Ax + 0,2) + (x + 0,2)? 

—(1+2x+x8]  (xER) 
=[x — 0,4x + 0,44] (xe€R) A 

Ejercicio 1 
Si fi, fe, fa, están definidas por 

f, :x——2x (xeR) 

fix——x? (xeR) 

Sai 5H (ER) 
hallar As(£,), Ao (f2) y Aso. (L5). Ed 
El operador A, aplica la función original tabulada f en Ar(f), y 
las imágenes de x,, x2,... bajo Ar(f) están tabuladas en la colum- 
na de las primeras diferencias. Por una segunda aplicación de A, 
podemos obtener la función A,(A,(f)); las imágenes de x;, x,... 


bajo esta función están tabuladas en la columna de las segundas di- 
ferencias. 


Nuestra forma diagramática anterior de la tabla de diferencias se 
trasforma ahora en 


1(x,) 


Ll ARAnt(x) 
Ah Aní(x) 
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Notación 3 


Ejemplo 1 


Ejercicio 1 
(3 minutos) 


MB4.2.1 


Usando la notación para composición de funciones introducida en la 
Sección 1.2.2 de la Unidad 1, Funciones, podemos escribir A, (A, 
(f)) = Are Anílf), que abreviaremos después en Añ(f) u, omitien- 
do completamente los paréntesis, en Af. Esta notación se puede 
extender de manera obvia para obtener las diferencias siguientes, de 
modo que, si denotamos la imagen del elemento x bajo la función A;f 
por Asfí(x), donde n E Z*, nuestra tabla general de diferencias to- 
ma ahora la forma 
t(x,) 


ps 
E Ant) 


N 


A 
rd 


"4 ES 
1(x,) ht) Abnt(a:) 
( a + 
O sá Anxo) Ani(a) 
e N 
160) CT 4160) 


N 7 
A di 
) 


(xs 


Las cantidades Asf(x) se llaman diferencias de orden n (para cual- Definición 3 


quier entero positivo n). El operador A; se llama operador diferen- .on$ + 
cia de orden n para la separación h. Definición 4 
" $* $ 
Ejercicio 2 Ejercicio 2 
(3 minutos) 


Si fi, fe y fa están definidas por 
Y 2x (xeR) 
Lun. (xe R) 
fix —= 3x2 +2x (xER) 


hallar AZ(J,), A3 (42), AUS). " 


Ejercicio 3 Ejercicio 3 
isis a w (5 minutos) 
Complete cada uno de los siguientes enunciados llenando las casillas 


como en (i): 


A,(x— senx) =| [x-——sen(x + h) — senx] 


(Todos los dominios no especificados son R o F según el caso.) 


6) ay —sema) > [A] 
(11) A,sen(x) = 
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El lector puede hallar difíciles los dos ejercicios siguientes: Si así su- 
cede, trate de entender las sugerencias y soluciones. No gaste mucho 
tiempo en ellos, pero en todo caso, observe el resultado de la parte (i) 
del primer ejercicio en la forma dada al final de la solución. 


Ejercicio 4 Ejercicio 4 
A , a , (10 minutos) 
(i) ¿Es A, compatible con la adición sobre F? Si su respuesta es 


SI, entonces vea si puede reconocer la operación binaria induci- 
da en el conjunto imagen. 

(ii) ¿Es A, compatible con la multiplicación sobre F? Si su respues- 
ta es SI, entonces vea si puede reconocer la operación binaria in- 
ducida en el conjunto imagen. 


(Para las definiciones de compatibilidad y de operación binaria indu- 
cida, véase la Unidad 3, Operaciones y morfismos). A continuación 
del ejercicio 5 se dan algunas sugerencias relativas a este ejercicio. M 


Ejercicio 5 Ejercicio 5 
(10 minutos) 


Si fy g E F son tales que Ar(f) = Ar(g), mostrar que 


a 
donde p es una función periódica con período h*. ¿Qué más se pue- 
de decir acerca de p si A,(f) = Ar(g) para todo h E R? (E 
Sugerencias para la solución del ejercicio 4 Sugerencias 


Para desarrollar un problema de este tipo que requiere un resultado 
muy general, es útil considerar algunos casos especiales simples 
(Polyat páginas 190-199). 


Por ejemplo, las siguientes dos funciones 


Fai (xeR) 


AA 1 (x€ R) 
tienen la misma imagen: 


AN): ——>2hx + h? (xeR) 


AMS): x—— 2hx + h? (x€R) 
como puede verificarse fácilmente. 


*Una función f con dominio R se dice que es una función periódica con período h Definición 5 

Fx + h) = f(x), para todo x E R. El gráfico de una tal función “se repite” a sí mis- * 

mo a intervalos de longitud h. Por ejemplo, la función seno es periódica con período 2r. 
+G. Polya, How to Solve it, Open University ed. (Doubleday Anchor Books 1970). Este 

libro es básico para el Curso Básico de Matemática; se hace referencia a él en los tex- 

tos como Polya. (continúa en la página 27) 
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Solución 1 
ASÍ) => Nx +5) f(x) (xeR) 
= [x—>Ax+5)-2x (xeR) 


= [x—— 10 (xe R)] 

Ao fa) = [lx —— L,4x + 0,49 (xe RJ] 
> 

A 0,01 (4) = [x—0 (xe RJ 


Solución 2 
A partir de la solución del ejercicio 1 tenemos 


As(f,) = [x—— 10 (xe RJ] 
Por tanto 


As(f1) = AsíAs/,) 

= [x—> 10 — 10 (xe RJ 
= [x-—>0 (xe R)] 

También a partir de allí tenemos 
Aor(f2) = [x——1,4x + 0,49 x€ R)] 

Por tanto 

Añ1(f2) = [lx —— (1,4(x +0,7) + 0,49; 
— (1,4x + 0,49) (x€ R)] 


= [x—— 0,98 (x € R)] 
Finalmente, 
Af) = [lx ——(3Ax + 1)? + Ax + 1) 
—(3x? + 2x) (x e R)] 
= [x——(3x? + 9x? + 9x + 3 + 2x +2) 
—(3x? + 2x] (x € R)] 


= [x—>9x? +9x+5 (xeR) 


Af) = [lx —> (Ax + 1)? + Ax + 1) + 5) 
—(9x? +9x +5) (xeR) 
= [x—>18x +18  (xeR) 
Solución 3 
(ii) sen (x + h) —senx 
(i1i) sen (2 + h) —sen2 
(iv) AM(x—— senx) = Ay(x——sen(x + h) —sen x) 
= [x——sen(x + 2h) — sen(x + h)] 


—[x— sen x + h) —senx] 


= [x——>sen(x + 2h) — 2sen(x + h) +senx] 


(v) AJu— x?) = Ax —— (x + h)? — x?) 


= Ax 3x?h + 3xh? + h?) 


= Ax — (3HMx + hy? + 3h (x + h) + h?) 


—(3x?h + 3xh? + h*)) 
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Solución 1 


Solución 2 


Solución 3 


= Ax 6h?x + 6h?) 
= [x—>(6h%(x + h) + 6h%) — (6h?x + 6h*)] 
= [x-— 6h?) [] 


(viene de la página 25) 


También 

CS (xe R) 
y 

82: >x +2 (xeR) 


tienen la misma imagen: 


Aj(g1):x=—h (xe R) 


A(g2):x—— h (xe R) 


De modo que, en lugar de considerar la compatibilidad en general, 
veamos primero si funciona para estos casos especiales, por ejemplo, 


¿es Af, +81) = Af + 82)? 


¿es ALÍ x 81) = Aí x 82)? 


Cuando haya respondido a estas preguntas, usted podría responder 
por lo menos a una de (i) y (ii). 


4.2.2 Diferencias de polinomios 


Cuando extrapoló la tabla de distancias de parada del “Código de 
carreteras” de la Sección 4.1.3, a 70 millas por hora, usó el hecho de 
que las segundas diferencias de esta tabla eran todas iguales. Cuan- 
do usted halló 37? extrapolando una tabla de cuadrados que termi- 
naba en 36?, usó de nuevo el hecho de que las segundas diferencias 
eran iguales. Quizás se pregunta si hay otras funciones para las cua- 
les las segundas diferencias sean todas iguales, o si las diferencias 
para una tabla de cubos tienen una propiedad similar. 


Siguiendo nuestro argumento semejante al del ejercicio 4.2.1.5, po- 
dríamos investigar funciones reales f, para las cuales Af sea una 
función constante para todo h E R, y tratar de deducir cómo es f. Es- 
te sería ciertamente un enfoque satisfactorio, pero podría ser un po- 
co difícil. Así que hemos escogido algo más fácil pero menos satisfac- 
torio y es pedir a usted que verifique el resultado que se da en el si- 
guiente ejercicio. 


Ejercicio 1 

Verifique que para la función 
f:ix—ax?+bx+c (xeR) 

(i) Arf es una función lineal, 


(ii) A?f es una función constante. E 
h 
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4.2.2 


Discusión 
* * 


Ejercicio 1 
(2 minutos) 


(continúa en la página 30) 


Solución 4.2.1.4 


Usamos la notación introducida en las sugerencias y tratamos prime- 
ro el caso (ii). 


AHÍ, x 21) = Ax —> x>) 
= [u—=>(x + MP 30) 
AHÍ x 82) = Ax —— x? + 2x? + x + 2) 
= [x—> ((x + hy + 2Ax +A) +(x + h) +2) 
=(x? + 2x? + x +2)) 


podemos ver fácilmente que Alf, X gi) 4 Anlfe X g2). De 
modo que este ejemplo muestra que A, no es compatible con la mul- 
tiplicación sobre F. 


Lo anterior ilustra un punto importante: siempre que queramos deter- 
minar si un enunciado general respecto de los elementos de un con- 
junto es verdadero o falso, entonces, si podemos hallar elementos 
cualesquiera del conjunto para los cuales el enunciado sea falso, en- 
tonces es falso. Por otra parte, si hallamos un elemento del conjunto 
para el cual el enunciado es verdadero, esto no prueba que sea verda- 
dero para todos los elementos del conjunto. En este caso nuestro enun- 
ciado general es 


“A, es compatible con la multiplicación sobre F” 


y este es un enunciado sobre todos los elementos de F. Hemos hallado 
cuatro elementos, fi, fe, £1, 82, para los cuales el enunciado es fal- 
so; por consiguiente, es falso. 


Por otra parte, como usted probablemente habrá hallado si siguió nues- 
tra sugerencia, el enunciado 


“A, es compatible con la adición en F” 


es verdadero para fi, fe, £1, £2. De modo que no hemos llegado a 
conclusión alguna sobre (i). Podemos continuar buscando elementos 
de F para los cuales sea falso el enunciado o tratar de probar que es 
verdadero, cualesquiera sean los elementos de F que escojamos. El ca- 
mino que adoptemos depende de nuestra intuición. ¿Creemos que es 
verdadero, o que es falso? Si nuestra intuición no está lo suficiente- 
mente desarrollada en este problema como para influenciar nuestro 
sentir, debemos ensayar otros ejemplos. 


Finalmente, en este caso, tendremos que sentarnos y probar que el 
enunciado es verdadero, como sigue: 


Sean fi, fe, 81, 82 E F tales que 
An) = An) 
Ar (21) = Ar(g2) 

Entonces 


AHÍ, +81) = Al f(x) + 8 ,(x)) 
= [pi fix +1) + 81(x + hi) 
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Solución 4.2.1.4 


Ecuación (1) 


Ecuación (2) 


(400 + 2,09) 
= [x—— f(x + h) — f(x) 
+ [x= g/(x +) - ¿09 
= Af) + Ay(8 1) 
De modo similar A,(f2 + g2) = Anlf2) + An(g2). Se sigue de las 
ecuaciones (1) y (2) que 


AHÍ, +21) = AL + 82) 


y como fi, f2, 81, $2 son elementos de F, se sigue que A, es com- 
patible con la adición sobre F. La operación binaria inducida [] sobre 
el conjunto imagen se define por 


Af) DA, (e) = Af + g) para todos los f, g E F 


y por nuestro trabajo anterior reconocemos [J como la adición de fun- 
ciones, 1.e. 


| AHÍ +8) = Af) + Ay(g) a 


Solución 4.2.1.5 
Ya que 
Af) = Ay(g) 
ANS) — Ajlg) = O 
donde O es la función constante 
x——>0 (xeR) 
Para una función real g, definimos —g como 
PR lo 
donde —g tiene el mismo dominio que g. Entonces podemos escribir 
ASÍ) — Ajlg) = AI) + (—Aj(g) = Af) + Ají—g) 
y por el resultado del ejercicio 4 tenemos 
AlSf+(-8)=0 o Alf-g)=0 
Ahora, si hacemos p = f — g, entonces A,(p) = O, así que tenemos 
píx + h) — píx) =0 


para todo x E R, de modo que p es periódica con período h, de lo cual 
se sigue el resultado requerido, a saber, f = g + p. 


Si la ecuación (3) es también verdadera para todo h E R, entonces, 
haciendo x = 0, tenemos 


píh) — p(0) =0 
así que 
plh) = p(0) 


para todo h E R,i.e. pes una función constante, cuya única imagen 
es p(0). (Para la definición de una función constante véase Unidad 1, 
Funciones.) E] 
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Solución 4.2.1.5 


Ecuación (3) 


(viene de la página 27) 


De este último ejercicio se sigue que una función general con la pro- 
piedad requerida de que sus segundas diferencias sean iguales es 


x—ax?+bx+c (xeR) 


No se sigue, sin embargo, que este sea el único tipo de función cuyas 
segundas diferencias sean constantes (cualquiera sea h), aunque es- 
to es en efecto cierto. Cualquier expresión de la forma ax? + bx +c 
(con a 4 0) se llama un polinomio cuadrático en x y una función que 
aplique x a un polinomio cuadrático en x se llama una función cuadrá- 
tica. 


Ejercicio 2 


El gráfico de una función cuadrática es una parábola. Si no está fa- 
miliarizado con esta curva, puede escoger valores numéricos para a, b, 
Cc y trazar una de tales curvas. | 


El problema que resolvimos se puede generalizar: ¿podemos hallar una 
función (real) general con terceras diferencias constantes o con dife- 
rencias de orden m constantes, donde m es cualquier entero positi- 
vo? Trate de adivinar la respuesta antes de proceder. Una conjetura 
plausible para una función general con terceras diferencias constantes 
es 


fix—oax +bx+cx+d (xeR) 


donde a, b, c, d son números reales cualesquiera. Cualquier expresión 
de la forma ax? + bx? + cx + d (con a 4 0) se llama un polinomio 
cúbicoen x y la función f es una función cúbica. Tomando las prime- 
ras diferencias de f tenemos . 


A,Í(x) = lalx + hy + b(x + hY +c(x +h) + d) 
(fax? + bx? + cx + d) 
= 3ahx? + (3ah? + 2bh)x + (ah? + bh? + ch) 


Es importante notar que el miembro de la derecha es un polinomio 
cuadrático en x; o sea que hemos mostrado que 


Ax (función cúbica) = (función cuadrática) 
Como ya sabemos, por el ejercicio 1, que 


Ar (función cuadrática) = (función lineal) 


A, (función lineal) = (función constante) 


se sigue que cualquier función cúbica tiene terceras diferencias cons- 
tantes. 


Este argumento se puede generalizar. La función (real) más general 
del mismo tipo que una función cúbica o cuadrática es 


¿9 0 ¡7 4 ¿+ 04x +4 (xER) 
y X 
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Definición 1 


* * 


Definición 2 


* $ 


Ejercicio 2 
(3 minutos) 


Tema principal 
- 


Definición 3 


* e 


Definición 4 


* * 


donde n es cualquier entero positivo 0 O y an, An-1,..., 41, y SON 
números reales cualesquiera, llamados coeficientes. 


Previsto que an *% 0, la expresión anx” + ::: + ay se llama 
un polinomio de grado n y la función correspondiente se llama función 
polinómica de grado n. Usted puede completar el argumento por su 
cuenta. 


Ejercicio 3 


¿Cuáles de los siguientes son polinomios? 


(1) x* — 6x +1, (11) (x — 3)M(x + 2), 


La ES . 

(111) ——, (iv) 2x +4, (v) 7. 
x-—3 

Otra vez, si encuentra difícil este ejercicio (o si no tiene tiempo) tra- 


te de entender la solución. En cualquier caso debe mirar el resultado 
del ejercicio. A 


Ejercicio 4 

Si f es una función polinómica de grado n, ¿es A,f una función po- 
linómica?, si lo es, ¿cuál es su grado? ¿Qué puede concluir acerca de 
la tabla de diferencias de una función polinómica, y en particular, 


acerca de la diferencia de orden n de una función polinómica de gra- 
do n? Si 


4.2.3 Resumen 


En esta sección hemos examinado el comportamiento de columnas su- 
cesivas de la tabla de diferencias para una función tabulada, y mos- 
tramos cómo pueden considerarse como generadas por aplicaciones 
sucesivas del operador de diferencia A», a la función original. Este 
operador proporciona una notación conveniente para expresar fórmu- 
las que contienen diferencias sucintamente. Se demostrará en una 
unidad posterior que tiene afinidades con el operador de diferencia- 
ción. A, es compatible con la operación de adición (y de sustracción), 
lo que es un resultado importante y útil. 


Después mostramos que las diferencias de orden n de una función po- 
linómica de grado n son todas iguales. Recíprocamente, si hallamos 
que las diferencias de orden n en una tabla de diferencias son todas 
iguales (cualquiera sea h), podemos suponer que la función que gene- 
ra la tabla es un polinomio de grado n. 
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Definición 5 
* * 


Definición 6 
h * $ 
Definición 7 
h »x $ 


Ejercicio 3 
(2 minutos) 


Ejercicio 4 
(2 minutos) 


4.2.3 


Resumen 
h * 


MB4.2.2 
Solución 4.2.2.1 Solución 4.2.21 


(1) ALÍ = [x——(a(x + hY? + b(x + h) +) 
—(ax? + bx + c) (xe R)] 
= [x—— 2ahx + ah? + bh (xeR)] 


que es lineal. 
(ii) ARS = [x—— (Qah(x + h) + ah? + bh) 
— (2ahx + ah? + bh) (xe R)] 
= [> 2ah* (xe R)] 


que es constante. E 


Solución 4.2.2.2 Solución 4.2.2.2 


Con b = c = 0, son ejemplos de curvas 


Para cualquier otra escogencia de a, b, c, se obtiene una curva simi- 
lar, desplazada del origen. a 
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Solución 4.2.2.3 


(i) es un polinomio (cuártico). 

(ii) es un polinomio (cuadrático). 
(iii) no es un polinomio. 
(iv) es un polinomio (lineal). 

(v) es un polinomio (constante). E 


Solución 4.2.2.4 
Fox +0 ¡7 4 ¿404x449 (xER) 
Aj ¿:x—— [aj (x + Y" +4, (x+h7*+---+a,(x + h) + a9) 
(a, + 4, ¡7 4---+4¡x +49) (xEeR) 
Ahora desarrollamos cada paréntesis (x + h)", k = 2,3,..., n, usan- 
do el desarrollo binomial, 


n n 
(x+ hy =x" + [oe + [ea Ho +h 


Así que Jan(x + h)” — anx” | es 
Lay + nx" 1h +---+h") — a,x") 
= a (nx"" 1h +--- + h") 
= polinomio de grado (n — 1) 


De manera similar podemos desarrollar cada uno de los términos y lue- 
go coleccionar todos los resultados. Pero lo que esencialmente debe 
observarse aquí es que comenzamos con f, una función polinómica de 
grado n y reconocemos a A,f como otra función polinómica, pero la 
potencia más alta de la variable x es ahora (n — 1), i.e. A,f es una 
función polinómica de grado (n — 1). 


Así que A, (función de grado n) = (función de grado n — 1). 


De modo similar 
A? (función de grado n) 


= A, (función de grado n — 1) 
= (función de grado n — 2) 
y 
A; (función de grado n) 
= (función de grado n — n) 
= (función constante) E 
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Solución 4.2.2.3 


Solución 4.2.2.4 


4.3 INTERPOLACION NO LINEAL 


4.3.0 Introducción 


En uno de los ejercicios que hizo antes en esta unidad (ejercicio 
4.1.2.3) vio que la interpolación lineal tendía a ser más exacta entre 
más pequeño fuera el espacio de separación entre los puntos tabula- 
res. Los diseñadores de buenas tablas de funciones usadas común- 
mente, tales como las de logaritmos y funciones trigonométricas, tie- 
nen cuidado de que la separación tabular sea suficientemente pequeña 
para que la interpolación lineal sea exacta y conveniente. Para una 
tabla que no se usa mucho, sin embargo, el diseñador puede muy bien 
decidirse a economizar calculando el valor de la función para relati- 
vamente pocos puntos tabulares; tiene particular tendencia a hacer 
esto si cada valor tabular demanda una cantidad de esfuerzo de cálcu- 
lo y si la tabla resultante se usa solamente una vez porque surge de 
una situación única, como en el caso de la función mencionada en la 
introducción a esta unidad que describe el efecto elevador del aire so- 
bre partes particulares del ala de un aeroplano bajo condiciones parti- 
culares de vuelo. En tales casos la separación tabular tiende a ser 
demasiado grande para justificar la interpolación lineal y se necesi- 
tan métodos más refinados. 


El propósito de esta parte de la unidad es mostrar al lector uno de ta- 
les métodos, en el cual la función tabulada se aproxima entre los pun- 
tos tabulares no por una función lineal sino por una función polinómi- 
ca. Usted se puede preguntar “¿Qué hay de especial en las funciones 
polinómicas? ¿Por qué no usar alguna otra función tal como seno o co- 
seno?” Una razón es simplemente que la interpolación lineal es muy 
conveniente, en gran parte porque las diferencias de polinomios tie- 
nen propiedades favorables, como se vio en la sección precedente; y 
la conveniencia es un factor muy importante al hacer una escogencia 


de métodos numéricos. 


Una razón más para usar polinomios está contenida en un teorema 
probado por un matemático alemán, Weierstrass, quien vivió en el si- 
glo pasado. Este teorema (cuya prueba no consideraremos aquí) ex- 
presa que cualquier función continua* puede aproximarse sobre un 
intervalo dado con cualquier aproximación deseada, por medio de un 
polinomio. Desafortunadamente, el teorema de Weierstrass tiene la 
característica molesta de tantos de los llamados teoremas de existen- 
cia en matemática: nos dice que el polinomio existe, pero no cómo ha- 
llarlo. En el resto de la sección no usaremos el teorema de Weierstrass 
excepto para confiar en la validez definitiva de los métodos de aproxi- 
mación polinómica, y nos concentraremos más bien en métodos prác- 
ticos para hallar polinomios aproximantes y para interpolar con su 
ayuda. 


*El término continua en este contexto significa que el gráfico de la función es conti- 
nuo: i.e. no tiene rompimientos o saltos la curva. El concepto se estudiará más rigu- 
rosamente en la Unidad 7, Sucesiones y límites 1. 
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4.3 


4.3.0 


Introducción 
* 


Karl W. T. Weierstrass -(1815-1897) 


4.3.1 Ajuste de una función lineal 


El tipo más simple de interpolación polinómica es el que hemos consi.- 
derado en la Sección 4.1.2; en efecto, las funciones lineales usadas 
allí para aproximar las funciones tabuladas son, por definición de 
“función polinómica” dada en la Sección 4.2.2 precisamente funcio- 
nes polinómicas de grado 1. Antes de discutir interpolación no lineal, 
por tanto, recapitulamos los principios de la interpolación lineal usan- 
do un punto de vista ligeramente diferente del que se usó en la sec- 
ción dedicada a este tópico. 


En la interpolación lineal aproximamos la función entre dos puntos 
tabulares sucesivos, digamos x» y xx+1, por una función lineal que 
Se ajusta 2 (i.e. tiene los mismos valores tabulares que) la función 
tabulada en esos dos puntos tabulares. (Véase la figura). 


Esta función lineal se puede denotar por /, donde 
lx ax +b (x€ [Xx Xz4 1)) 


y a, b son dos números reales. 


Ejercicio 1 


¿Por qué usamos [xx, x;+,], en lugar de R, para el dominio de 1? 
E 


Hemos estipulado que | debe ajustarse a la función tabulada en xx y 
Xx+13 O sea 


xy = f(x ) 


Uxr+1) = FXe+ 1) 


Estas dos condiciones son suficientes para determinar las dos cons- 
tantes a y b de la definición de l, pues son equivalentes a 


xa +b=f(x,) 
Xx+19 + b= f(Xi41) 


y como los puntos tabulares xs, x1+1 y los valores fíxx), fíxx+1) 
se conocen, podemos tratar estas ecuaciones como una pareja de ecua- 
ciones simultáneas para las dos incógnitas a y b. 
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4.3.1 


Discusión 
* * 


Definición 1 
k + 


Ejercicio 1 
(2 minutos) 


(continúa en la página 36) 


Solución 1 
El gráfico de l va a ser un segmento de recta: no una recta completa. 
El equivalente algebraico de esta condición geométrica es la restric- 


ción sobre el dominio de /. Evidentemente queremos aproximar a la 
curva completa como se ve a continuación. 


no como esto: 


(viene de la página 35) 


Ejercicio 2 


Resolver las ecuaciones para a y b y obtener una fórmula explícita pa- 
ra l(x) en términos de x y de los puntos y valores tabulares. 


(La solución de este ejercicio requiere cierto grado de facilidad en la 
manipulación algebraica. Si no tiene esta facilidad, lea la solución y 
luego trabaje por sí mismo en los detalles, paso por paso, si tiene tiem- 
po. Lo importante es tomar nota del resultado de la ecuación (3). En 
general, en esta sección, los detalles técnicos (algunos de los cuales se 
omiten en cualquier caso) no son tan importantes como su compren- 
sión de lo que estamos tratando de hacer y una apreciación de los 


resultados.) En 

Ejercicio 3 

Use la ecuación (3), página 38, para hacer un es- Xx tan x 

timativo de tan (1,444) a partir de la tabla de la 

derecha, y comparar con el verdadero valor (con 1,43 7,055 

tres cifras decimales), tan (1,444) = 7,844. 1,44 7,602 
1,45 8,238 
1,46 8,989 


(Usted puede, desde luego, responder a este ejercicio por los métodos 
de la Sección 4.1.2, pero, para lograr familiarización con los métodos de 
ahora, le sugerimos usar la ecuación (3).) pal 
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Solución 1 


Ejercicio 2 
(3 minutos) 


Ejercicio 3 
(3 minutos) 


4.3.2 Interpolación polinómica 


Cuando la interpolación lineal no es muy exacta, como en el ejercicio 
anterior, ello indica que la función tabulada no está muy exactamen- 
te representada por la función lineal /(x) sobre el intervalo [xx, xx+1). 
En tales casos podemos dar paso a la no linealidad usando una fun- 
ción aproximante simple no lineal en lugar de una lineal. Si la función 
aproximante va a ser una polinómica, entonces la más simple (i.e. de 
menor grado) posibilidad no lineal es la función cuadrática 


q:ix—axd+bx+c 


Para completar la especificación de esta función necesitamos valores 
para las tres constantes a, b, c y también una especificación del domi- 
nio. Siguiendo el método usado en la sección anterior podríamos tra- 
tar de determinar a, b, y c por la condición de que q se ajuste a la 
función tabulada en dos puntos tabulares sucesivos, digamos x» y 
xx+1 ¿Puede ver usted la dificultad que surge? Obtendríamos dos 
ecuaciones simultáneas 


xa + xb+c= f(x;) ] 


Xi 10 + Xp 1D + 0 = f(Xy4 1) 


pero no son suficientes para determinar a, b, y c (véase figura). 


Para hallar tres cantidades desconocidas necesitamos tres ecuaciones 
simultáneas. Hay varias maneras de obtener una tercera ecuación; la 
más simple y la que consideraremos aquí, es imponer a q la condición 
de que se ajuste a la función tabulada en tres puntos tabulares con- 
secutivos en lugar de solo dos (véase figura). 


Gráfico de la única 
cuadrática que se ajusta a 
la función tabulada en 

los 3 puntos tabulares 
consecutivos. 


Gráfico de una de las 
muchas cuadráticas 

que se ajustan a la 
función tabulada 

en solo 2 puntos tabulares. 


flx..2) 
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4.3.2 


Discusión 
Rh * 


(continúa en la página 39) 


MB4.3.1 
Solución 4.3.1.2 Solución 4.3.1.2 


La solución de las ecuaciones simultáneas es 


eq Fr 1) — Sa) 
Xk+1 — Xx 
De Xu+1 (Xp) — Au (Xi 1) 
Xx+a TX 
que sustituida en la fórmula de l(x) da 


Xx AS o e 
E Fx) + —— JXx+ 1) Ecuación (1) 
Xk+1 — Xx Xx+ 1 — Xx 


1) = 


Trasformaremos esta solución un poco más y eliminaremos x en favor 
de la variable 


que representa la proporción en la cual el número x divide al intervalo 
[x», xx+1] (véase figura). 


o :] 1 
—_—_—> __—_—_. >= >_—————. 
¡a x Y Xues 
«———— h —_—s 
Observando que 

Xiri —= X Xx 

li 3 | = 10 

Xx+1 7 Xx Xu TX 


la ecuación (1) se trasforma ahora en 


lo) = f(x) + A 1) — FOXIJO, (E [Xi Xi + 1)) Ecuación (2) 


Una manera más elegante de escribir la ecuación (2) es usar el opera- 
dor diferencia definido en la Sección 4.2.1; se simplifica entonces en 


Ux) = f(x) +0.A, f(x) (xElXi Xx+ 1]) WM Ecuación (3) 
Solución 4.3.1.3 Solución 4.3.1.3 
A la derecha se muestra la tabla de di- 
ferencias. El intervalo [xx, xx+1] que AS tanx  A,p (tan x) ] 
contiene 1,444 es [1,44, 1,45]> por con-. === 
siguiente xz = 1,44 y 1,43 7,055 0.547 

1,44 7,602 0.636 
0,004 4 1,45 8,238 0.751 
= 01 710 1,46 8,989 


La ecuación (3) da entonces 

tan (1,444) = 7,602 + $5 X 0,636 = 7,856 
que es el estimativo de tan (1,444) dado por interpolación lineal. El 
error de 0,012 muestra que la aproximación de tan x por l(x) no da bue- 


na exactitud en este intervalo. Bm 
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(viene de la página 37) 


Llamando a estos puntos Xx», Xx+1, Xx+2, tenemos entonces las tres 
ecuaciones 


q(x,) = xía + xb + <= f(x,) 
a: PS 
Aur 1) = Xi 1 + Xp 1D + 0 = f(Xu4 1) 
2 
AXx+2) = Xi 20 + Xp+2b + 0 = f(Xp4 2) 
y como el número de ecuaciones es ahora igual al número de incógni- 
tas, las ecuaciones contienen ahora suficiente información para deter- 
minara, b yc. 
Queda por especificar el dominio de q(x). Por las ecuaciones (1) el 
dominio debe incluir por lo menos los números Xx, Xk+1) Y Xk+2, Y 
si la función debe utilizarse para interpolación, debe también incluir 


los valores intermedios; así que el dominio mínimo para q(x) es [xx, 
Xx+2]. El gráfico de esta función 


qx —ax? +bx+c (x€[Xi Xx+2)) 


es el que se muestra en la figura anterior. Desde luego, si se necesita- 
ra la función (y fuera adecuada) para extrapolación, el dominio se po- 
dría extender con facilidad. 


Para hallar las constantes a, b, c, el método obvio de proceder es re- 
solver las tres ecuaciones simultáneas (1). 


Hay técnicas establecidas para resolver tal sistema de ecuaciones; 
usted hallará algunos de ellos en una de las unidades sobre Algebra 
lineal. Aquí, sin embargo, no tenemos que ver con las técnicas sino so- 
lamente con el resultado final, que puede ponerse en la forma 


qx) =ax?+bx+c 
Am XD — Mg 4 2) 
(a — XX a JA 


(e = AJA = Fxa:2) 
e => y A+ 


(x — xx — Xg4 1) 
(Xx+2 — XX +2 — Xx +1) 


+ FXs+ 2) 


La expresión de la derecha se llama interpolación polinómica de La- 
grange; es la análoga de la ecuación (1) de 4.3.1 para el caso lineal. (No 
es necesario que se aprenda esta fórmula de memoria.) 


Como en nuestro tratamiento del caso lineal, podemos simplificar la 
ecuación (2) usando la variable 0, relacionada con x como en la fi- 
gura: 


Entonces la ecuación (2) se trasforma (después de algunos detalles 
técnicos que se omiten) en un polinomio cuadrático en 6: 


qx) = 0 f(Xu+ 2) — 24 (X1 41) + SO) 
+0 —4/(X1+2) + 28141) — 3/0) + $ (Xu) 
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MB 4.3.2 


Ecuaciones (1) 


Joseph-Louis Lagrange (1736-1813) 


Ecuación (2) 


Definición 1 
* $ 


Ecuación (3) 


Esta fórmula no tiene mucha importancia; se escribe solamente para 
que usted pueda ver si reconoce el coeficiente 0? (i.e. el número que 
multiplica a 9? en el primer término de la derecha). Aparte del factor 
jeste coeficiente es Affí(xx). El coeficiente de 0 parece un poco 
más complicado pero puede escribirse en términos del operador dife- 
rencia; es Arflxr) — 3A5f(x1). Usando estas expresiones para los 
coeficientes podemos escribir la ecuación (3) en la forma concisa 


a) = f(x) + OA, f(x,) + 3010 — DAS (1) (xEl[Xz> Xu +2)) 


Este es un caso especial de una fórmula general llamada Fórmula de 


Gregory-Newton. 


Ejemplo 1 
Como ilustración de cómo usar la fórmula de Gregory-Newton, consi- 
deremos de nuevo la tabla del “Código de carreteras”, ¡.e. 


E E y A a ls MA 
Distancia de 


Rapidez parada Primeras Segundas 
(millas/h) (pies) diferencias diferencias 
20 40 
30 75 pa 10 
40 120 ss 10 
S0 175 6s 10 
60 240 


(«EKG oo 


Estamos interesados en un estimativo de la distancia en pies cuando 
se viaja a 24 millas/h. 


En este caso deberíamos escoger x» = 20 y Xx+1 = 30. Entonces 
g= 7% ¿A 
XxX: 30-20 10 
y tenemos h = 10; La ecuación (4) se trasforma en 
(24) = f(20) + A XA ¡0f(Q0) + - x e x mo. A?,f(20) 
...o > 210 To 
140 120 
"+50 
= 40 + 14— 1,2= 52,8 pies le] 


Ejercicio 1 


Usar la fórmula de interpolación cuadrática de 


Gregory-Newton para calcular tan (1,444) a partir x tan x 
de la tabla de la derecha y comparar con el resul- 
tado obtenido por interpolación lineal en el ejer- 1,43 7,055 
cicio 4.3.1.3 y con el resultado correcto 7,844. 1,44 7,602 
1,45 8,238 
1,46 8,989 
E 
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Ecuación (4) 


Definición 2 
€ *€* + 


Ejemplo 1 


Ejercicio 1 
(5 minutos) 


4.3.3 Fórmulas a partir de tablas 


Aunque la principal aplicación de la fórmula de Gregory-Newton es a 
la interpolación no lineal, puede usarse también para lo que esencial- 
mente es extrapolación no lineal de funciones polinómicas. Si el pro- 
pósito de una extrapolación es solamente hallar el elemento siguiente 
de una sucesión de números tales como 1,3, 5,7,... 640, 75, 120, 175, 
240,... (las distancias de parada del “Código de carreteras” conside- 
rado antes en esta unidad), entonces no hay necesidad de la fórmula 
de Gregory-Newton; simplemente construimos la función directamen- 
te a partir de la tabla de diferencias como se explicó en la Sección 
4.1.3. La fórmula de Gregory-Newton se vuelve útil, sin embargo, 
cuando queremos una fórmula que dé todos los elementos no tabula- 
dos de la sucesión. Por ejemplo, podemos usar la fórmula de Gregory - 
Newton para obtener una fórmula general que dé el elemento de orden 
n de la sucesión cuyos primeros cinco elementos son 


121242 124243,124+243 4421242? 4 3? 
A 


o sea, que nos dé una fórmula para la suma de los cuadrados de los 
primeros n números naturales. Tales fórmulas para suma de potencias 
se usarán en la Unidad 9, Integración 1 para calcular áreas. Un caso 
típico al cual se aplican los resultados es el que se ilustra en la figura; 
el área entre la curva y el eje x es aproximadamente la de un conjunto 
de rectángulos. 
y 
Gráfico de y=x(1-x) 


O 0,5 1 x 


La anchura de cada rectángulo es +5; la altura del rectángulo que ocu- 
n n 
10 l- al de modo que su 


La suma de las áreas de los 10 rectángulos se pue- 


pa el puesto n a partir de la izquierda es 

n n? 
100 1000" 
de calcular usando fórmulas deducidas en esta sección. (Para más de- 
talles véase la Unidad 9, Integración 1.) 


área es 


Para usar el método, consideremos otra vez la tabla de distancias de 
parada. 


S d Diferencias 
Rapidez Dist. de parada : 
(mi/h) (pies) Primera Segunda 

20 40 

30 75 E 10 

40 120 55 10 

50 175 65 10 

60 240 
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4.3.3 


Discusión 
h + 


(continúa en la página 42) 


Solución 4.3.2.1 


A la derecha se muestra la tabla de 


diferencias. El número 1,444 cae en De tan x A A? 

el intervalo [x», xx+2] con xx. = 1,43 

Ó 1,44. Se puede usar cualquier va- 1,43 — 7,055 547 

lor de x4. 1,44 7,602 636 89 
1,145 8,238 751 115 
1,46 8,989 


e 


Xx = 1,43 X; = 1,44 
0'=1,4 6 = 0£ 
q(x) = 7,055 + 1,4 x 0,547 q(x) = 7,602 + 0,4 x 0,636 
1,4 z 0,4 x 0,089 interpolación lineal 
¿04 pS M8, 0,115 
2 
=7,055 + 0,766 + 0,025 = 7,602 + 0,254 + (—0,014) 
= 7,846 = 7,842 | 
Escogiendo x+ = 1,43 ó 1,44, nuestra solución está de acuerdo con el 


valor correcto a menos de 0,002. Usando interpolación lineal (véase 
página 38) hemos hallado que tan (1,444) = 7,856, que no es muy exacto. 


La fórmula de Gregory-Newton en forma general para interpolación 
polinómica de grado n a la función f es 


g(%) = (xi) + 08, f(x1) + $0(0 — MARS (xy) + 
+00 = 1)(0 — 2)---(0—n + 1A5f(x) 


El dominio sobre el cual es útil para interpolación es el intervalo 
[xx, xx+n]. Es la fórmula básica para interpolación polinómica 
(aunque no la más práctica cuando el grado del polinomio es superior 
a 2). 


o 


(viene de la página 41) 


¿Podemos deducir a partir de la tabla una fórmula polinómica que dé 
la distancia de parada d en términos de la rapidez s? Esto es, ¿pode- 
mos hallar una función polinómica p tal que 


d = pís) 


sea válido para todos los puntos tabulares? 


Las primeras diferencias no son constantes; por consiguiente p no es 
de grado 1. Por otra parte, las segundas diferencias tabuladas son 
constantes, de modo que los valores tabulares se pueden ajustar por 
un polinomio de grado 2, i.e. cuadrático. La tabla no prueba que el gra- 
do del polinomio sea 2, ni que la función s——d es un polinomio. Pe- 
ro existe un polinomio cuadrático que se ajusta a todos los puntos ta- 
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MB4.3.2/4.3.3 
Solución 4.3.2.1 


Definición 3 
h »* $ 


bulados y la tabla proviene de una ley física, de modo que es un buen 
método científico adoptar el cuadrático como hipótesis de trabajo has- 
ta que aparezca una buena razón para sustituirlo (tal razón existiría, 
por ejemplo, si uno fuera lo suficientemente imprudente como para 
tratar de extrapolar a velocidades supersónicas o aún a velocidades 
de 150 mi/h). Véase también la discusión de la Sección 4.1.1. 


Una manera de hallar la función polinómica p sería sustituir tres va- 
lores tabulados diferentes de s en la fórmula 
d = p(s) = as? + bs+c< 


junto con los valores correspondientes de d y resolver las tres ecua- 
ciones simultáneas resultantes para a, b, c. Es mucho más rápido, sin 
embargo, usar el resultado de la Sección 4.3.2, donde'se resolvió un 
problema general de este tipo y nos dio la fórmula (4) de 4.3.2 de Gre- 
gory-Newton. En la notación actual esta fórmula es 


0(0 — 1 
dis) = dí5) + GA10 dis) + 5 Modo) (sEl5, +2) 
donde 0% está definida por 


s-S 
0 = z 


Ski+1 — Sk 
No importa a qué punto tabular llamemos sz, mientras conozcamos 
dis»), Ajo dís:) y Atodísr*). Tomemos sz = 20. Con esto y los 


datos de la tabla, la ecuación (1) se trasforma en 


d = 40 + 350 + pt (s e [20, 40]) 


= 40 + 300 + 50? 


con 


Sustituyendo 4 a partir de la ecuación (3) en la ecuación (2) tenemos 


E 208 
S 8 e 20) 


d =40 + 307 e 


2 
S 

4 =60+- 

0 + 3s + 2s + 20 


y? 
s + 35 (se[20,40) 


La función polinómica s——>d definida por la ecuación (4) se ha esco- 
gido para ajustar los valores tabulares sobre el intervalo [20, 40], 
y como solamente un cuadrático se ajusta a los tres puntos tabulares, 
es la única función cuadrática que lo hace. En consecuencia, si que- 
remos ajustar a los valores tabulares un cuadrático sobre el intervalo 
tabular, o uno mayor, debemos usar de nuevo el polinomio de la ecua- 
ción (4). Esta es entonces nuestra fórmula de extrapolación: 
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Ecuación (1) 


Ecuación (2) 


Ecuación (3) 


Ecuación (4) 


los puntos a través de los 
cuales pasa la curva porque 
las segundas diferencias son 
constantes 


240 


200 


160 


120 


los puntos a través de 
los cuales pasa la curva 
por construcción 


Gráfico de d 2... 


80 


eje s 


Sabemos, por la manera en que fue calculado, que el polinomio de la 
ecuación (5) se ajusta a los valores tabulares en el intervalo [20, 40], 
i.e. para s = 20, 30 y 40. Por la constancia de la columna de las segun- 
das diferencias sabemos que hay un cuadrático que se ajusta a todos 
los valores tabulares y por la unicidad de la fórmula de Gregory- New- 
ton, que este cuadrático debe ser el de la ecuación (5), de modo que 
este polinomio se ajusta a los valores tabulares s = 50 y 60 también. 
Para valores de s distintos de esos cinco la tabla no nos dice si el po- 
linomio de la ecuación (5) es una fórmula exacta para la distancia de 
parada; pero, como se explicó antes en esta sección, la fórmula es una 
buena hipótesis para confrontar con cualquier información ulterior 
que podamos tener sobre la fuente de los números de la tabla. Como 
ocurre en el caso presente, la “información ulterior” hasta cierto pun- 
to respalda la hipótesis. Si observa de nuevo la tabla de distancias 
de parada, verá que las distancias de parada dadas se clasifican en 
“distancias pensadas” y “distancias de frenado” (véase tabla si- 
guiente), 


Distancia 
Rapidez 
Parada Pensada Frenado 
20 40 20 20 
30 TS 30 45 
40 120 40 80 
S0 175 50 125 
60 240 60 180 


y estas dos distancias corresponden exactamente a los dos términos 
del polinomio de la ecuación (5). Como última parte del trabajo de de- 
tective, usted puede determinar por su cuenta qué hipótesis físicas se 
podrían adoptar para justificar los números particulares usados por el 
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Ecuación (5) 


ministro del trasporte para las dos componentes de la distancia de 
parada. 


Volvamos al problema mencionado antes: hallar una fórmula para la 
suma de los cuadrados de los primeros n números naturales. El mé- 
todo de resolverlo es idéntico al del ejemplo de las distancias de para- 
da, de modo que podrá hacerlo usted mismo con poca ayuda. 


Ejercicio 1 


Si se representa con S2(n) la suma de los cuadrados de los n prime- 


ros números naturales, ¡.e. para 1? + 2? + --- + n?, tabular el va- 
lor de S2(n) para n = 1, 2, 3, 4, 5, 6. ¿Cuál es el grado mínimo de un 
polinomio en n que se ajuste a esos valores? E] 


Ejercicio 2 


Hallar el polinomio de grado mínimo que se ajuste a S2(n) para los 
valores tabulados en su solución anterior. Verificarlo considerando 
S2(7). (En lo que sigue demostraremos que la fórmula para S2(n) 
es correcta para todo entero positivo n.) ES 


Ejercicio 3 
Hallar 


(i) la suma de los primeros n números naturales y 
(ii) la suma de los cubos de los primeros n números naturales. ll 


(Si está corto de tiempo, no gaste mucho en este ejercicio; léalo y en- 
tienda la solución y observe los resultados.) 
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Ejercicio 1 
(3 minutos) 


Ejercicio 2 
(10 minutos) 


Ejercicio 3 
(10 minutos) 


Solución 1 


n S2(n) A A? A? 
1 1 
2 5 E 5 , 
E 14 7 

16 2 
4 30 25 9 ) 
5 9 36 11 
6 91 


Como la columna de las terceras diferencias es la primera que es cons- 
tante, el polinomio es por lo menos del tercer grado. | 


Solución 2 
La fórmula de Gregory-Newton con x» = 1 queda: 


sim =1+40+ 5%, 290-0072) 


conó = n — 1. 
Simplificando se obtiene: 


6 + 24(n — 1) + 15(n — 1)(n — 2) + 2(n — 1)(n — 2)(n — 3) 
Sn)  —_———>>=———2 2 


6 
_ 6 +24n — 24 + 15n? — 45n + 30 + 2n? — 12n? + 22n- 12 
> 6 
_n+3n* +2n* 
sa 6 
o sea 
n(n + 1)(Qn + 1) 
S2(n) = A ape a 
Solución 3 
(i) Sim=1+2+3+---+n 
La tabla de diferencias para S,(n) es 
n S,(n) A¡Si(n)  A¿S,(n) 
1 1 
2 3 : 1 
3 6 4 1 
4 10 5 1 
S 15 6 1 
6 21 


Como las segundas diferencias son constantes, el polinomio es de 
segundo grado (por lo menos). 


La fórmula de Gregory-Newton para el polinomio cuadrático Sí(n) es 
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Solución 1 


Solución 2 


Solución 3 


Sin) = S(n,) + 9A,S(n,) + 30(0 — DATS(n) 


Tomamos nz = 1y 0 = ———=>—=n-1 


Se sigue que 
Si) =1+(M1—=1)x24+ Kn — Dn —2) Xx 1 
= inín + 1) 
(ii) Sinm=10+24+34+--40 


La tabla de diferencias para S3(n) es: 


n Sin)  A¡Sín) A?Sn) AjSam) AiSyn) 
1 1 
8 
2 9 19 
3 36 na 37 ln 6 
4 100 ¡3 61 sá 6 
3 225 216 91 
6 441 


Sa(n) = Sín,) + 9A,Sí(n,) + $00 — DAÍS(n,) 
+10(0 — 1)(0 — 2)A¡S(n,) 
+240(0 — 10 — 2X0 — 3)A1Sx(n,) 
Tomamos de nuevo nz =1y0 = nm — ly así, 
Sim) =1+(n— 1) x 8 + Hn — 1)(n — 2) x 19 
+ Un — 1)(n — 2)(n — 3) x 18 
+2 (n — 1Xn — 2Xn — 3Xn — 4)x 6 


Esto conlleva una buena cantidad de cálculos que dan el siguiente 


resultado 


n? e 
Sn) = 30 + 1) 


Observación 
Hay un “truco” para simplificar los cálculos; poner 


Sin) =0 +19 +... +nm 


i.e. comenzar con n; = 0 y se tiene entonces 
0 — “S Es 
Sin) =0 + 19 +7% a ¡990 — DO - 2 
2 6 
¡PO HA con 0 = n 
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En su solución al ejercicio 2, observará que las primeras diferencias 
de S2(n) son los cuadrados de los números naturales y la función 
“elevar al cuadrado” es una función polinómica cuadrática; por con- 
siguiente S2(n) debe ser un polinomio de un grado mayor, o sea cú- 
bico. Su solución mostraba también que el único polinomio cúbice que 
se ajusta a los valores de S, (n) para n = 1, 2, 3, 4 es 


n(n + 1)(Qn + 1) 
6 


Para ser precisos deberíamos argumentar de la manera siguiente: 


Tenemos 


Sín+1)=12+22+---+ 1 +(n + 1) 


Sán) =12 +22 +---+n?; 
se sigue que 
A,S2(n) = San + 1) — Sym) ' 
= (n + 1)?, de modo que 
S2(n) = C(n) + pín) 


donde C(n) = polinomio cúbico en n y p(n) = función periódica de 
período 1, i.e.* 
p(1) = p(2) =--- = p(n) (Véase el ejercicio 4.2.1.5) 


Si hacemos ahora C(n) AE 


p(1), pero S2(1) = C(1) = 1, i.e. p(1) = 0 por tanto pín) = 0. De modo 
que finalmente 


, entonces S2(1) = C(1) + 


nn + 1)(2n + 1) 


S2(n) = C(n) = 6 


En consecuencia Sz(n) debe ser igual a este polinomio no sólo para 
n = 1, 2, 3 y 4 sino para todos los enteros positivos, y concluimos que 


_ hín + 1)(Qn + 1) 


PIP Pe 6 
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MB4.3.4 


4.3.4 Resumen 4.3.4 
En esta parte hemos mostrado que los valores interpolados y extrapo- Resumen 


* $ 


lados de una función tabulada se pueden calcular por la fórmula de 
Gregory - Newton 


0-1 
80) = flop) + 08,S050 + 430%) 
9 — 1)---(0 — 1 
a 5) 
donde 
g= 7% 
Xurr T Xx 


(con la previsión usual acerca del peligro de la extrapolación). 


La fórmula es la base de muchas que se usan mucho en cálculos prác- 
ticos, especialmente relacionados con un computador. (Sin embargo, 
no es particularmente adecuada por cuanto con fórmulas alternas se 
puede obtener el mismo grado de exactitud con una cantidad de cálcu- 
los considerablemente menor.) Un segundo uso de la fórmula es para 
deducir fórmulas para la suma de los elementos de una sucesión de 
números, tal como la suma de potencias de los números naturales. Es- 
tas fórmulas se usarán en una unidad posterior que trata de integra- 
ción y ocasionalmente se necesitan en otros campos. 
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4.4 ERRORES 


4.4.0 Introducción 


Esta sección da otra aplicación a las diferencias finitas. Creemos que 
es suficientemente importante e interesante para incluirla aquí, pero 
como las unidades subsiguientes no dependen de esta sección, puede 
omitirla si está corto de tiempo. 


Hasta ahora hemos tratado los valores tabulares como valores exactos 
de la función en los puntos tabulares. Esta es, sin embargo, una super- 
simplificación. Los valores tabulares son casi siempre inexactos, por 
razón nada más que de los errores de redondeo; además, en demasia- 
dos casos, los valores tabulares ocasionales no son ni aproximada- 
mente correctos por alguna equivocación que proviene bien de hallar 
el valor de las imágenes de la función o de copiar los números en la 
tabla. En esta sección aprenderá cómo reconocer y tener en cuenta el 
efecto de errores de redondeo cuando se usa una tabla de diferencias, 
y cómo la tabla de diferencias hace posible, en casos favorables, loca- 
lizar y aún corregir una entrada tabular que es incorrecta por causa 
de una equivocación. 


4.4.1 Equivocaciones 


Consideramos las equivocaciones, primeramente porque sus efectos 
sobre la tabla de diferencias son más considerables que los de los erro- 
res de redondeo. He aquí dos tablas de valores de seis cifras para los 
cosenos de los números comprendidos entre 0,512 y 0,518 con una se- 
paración de 0,001. Una de las tablas contiene un número deliberada- 
mente errado, tal como podría ocurrir al copiar. Comparando las dos 
tablas puede verse cuál punto tabular tiene el valor tabular incorrec- 
to en una de las dos tablas; ¿pero, podría decir cuál tabla es la inco- 
rrecta? 


Tabla 1 Tabla II 
Xx cos x b ?% cos x 
0,512 0,871 766 0,512 0,871 766 
0,513 0,871 276 0,513 0,871 276 
0,514 0,870 785 0,514 0,870 785 
0,515 0,870 293 0,515 0,870 239 
0,516 0,869 800 0,516 0,869 800 
0017 0,869 306 0,517 0,869 306 
0,518 0,868 811 0,518 0,868 811 


Una manera de hallar si la tabla lo la Il es la que contiene el error en 
x = 0,515 es construir sus tablas de diferencias (la coma decimal no 
se tiene en cuenta usualmente al hacer esto): 
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4.4 


4.4.0 


Introducción 
* 


4.4.1 


Discusión 
* 


Tablas 1, II 


Tabla I con Diferencias 
Diferencias 
x ? cos x 


Primera Segunda 


0,512 0,871 766 


0,513  0,871276 e 2Í 
0,514  0,870785 400 =1 
0,515  0,870293 ¿o =1 
0,516 0,869800 ¿oy A: 
0,517 0,869 306 “i 


0,518 0,868 811 ld 


Tabla II con Diferencias 
Diferencias 
Xx 2 cos x 


Primera Segunda Tercera 


0,512  0,871766 


0,513  0,871276 Lo a 
0,514 0,870785 Tia, 55 ¿yg 
0,515 0,870239 Zas9  +107 0 
0,516 0,869800 gay 55 “o 
0,517 0,869306  4gs -1 


0,518 0,868 811 


En la tabla llas segundas diferencias son constantes, de modo que 
los valores tabulares se pueden ajustar por medio de un polinomio cua- 
drático en x; pero en la tabla II las últimas columnas de diferencias 
son muy irregulares, de modo que ningún polinomio simple da un ajus- 
te razonable a los valores tabulares dados. Como el gráfico de la 
función coseno es una curva más bien suave, esta es una evidencia 
muy fuerte de que la tabla Il es la que tiene el error y de que el valor 
0,870 239 dado en la tabla para x = 0,515 debe ser realmente 0,870 293 
como en la tabla L 


No solamente podemos usar las diferencias para determinar si hay o 
no un valor errado en la tabla, sino que podemos usarlas también para 
localizar el valor incorrecto y aun para corregirlo. El método es una ex- 
tensión de un método que los científicos experimentales usan. Al ha- 
cer una serie de medidas de alguna cantidad en una sucesión de 
instantes diferentes o de diferentes valores de alguna condición ex- 
perimental como la temperatura, es práctica común colocar los resul- 
tados en un gráfico. Si resulta que todos los puntos excepto uno caen 
sobre una curva suave, el experimentador toma esto como evidencia 
de que hubo un error al tomar la correspondiente medida y si es posi- 
ble repite la medición para ver qué estuvo mal. Si no se puede repetir 
la medida hay una fuerte tendencia a rechazar la medida anómala 
completamente, pero lo que se requiere realmente es alguna explica- 
ción de por qué es anómala; puede no ser un error, sino que puede in- 
dicar algún resultado físico significativo!) 
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Para ver cómo puede usarse la tabla de diferencias para detectar va- 
lores incorrectos o anómalos en la tabla, supongamos que una equivo- 
cación remplaza el valor que debería ser f(x»), en la tabla de f, por un 
valor erróneo que denotamos por fíxx) + E; así que E es el error in- 
troducido por la equivocación en el valor del punto tabular x, (véase 
la tabla ID. 


Xx f 


Xk-2 Fxi-2) 


Xk-1 Fxi-1) 
Xk Fx) + E 
Xk+1 Fis 1) 


Xk+2 FXi+ 2) 


Podemos considerar a ésta como la tabla de una función que es la su- 
ma de f y una “función error” cuyo valor tabular en xx es E y en to- 
dos los otros puntos tabulares es cero. En vista de que A; es un ho- 
momorfismo con respecto a la adición, ¡.e. que 


Arílf + función error) = A,f + Ar (función error) 


la tabla de diferencias asociada con la tabla III será la “suma” de la 
tabla de diferencias de f y la de la “función error” que es 


Esto muestra que el error E afecta más y más entradas en cada co- 
lumna sucesiva de diferencias; las entradas afectadas forman un es- 
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Tabla HI 


Tabla IV 


quema rectangular. En cualquier columna particular, el mayor efecto 
ocurre en las entradas cercanas a la entrada original errónea. (Los va- 
lores de los errores introducidos en la r-ésima columna de las di- 
ferencias, son proporcionales a los coeficientes del desarrollo bino- 
mial de (a — b)”.) 


Para ilustrar el método de localizar y corregir el valor incorrecto apli- 
quémoslo a la tabla II, que está repetida a continuación con sus co- 
lumnas de diferencias. 


Tabla II con Diferencias 


Diferencias 
Xx 2 cos x 
Primera Segunda Tercera 


0,512 0,871 766 490 


0,513  0,871276 -1 
0,514  0,870785 Mes -55 pu 
0,515 0,870239 “yzg +107 T¡¿ 
0,516 0,869800 _y9y 55 sg 
0,517 0,869306 "¿os -1 


0,518 0,868 811 


En la columna de primeras diferencias, es difícil desenmarañar la 
contribución del error E a partir de la contribución de la función ta- 
bulada, pero en la segunda y tercera columnas la contribución del error 
E predomina y se puede discernir claramente el esquema que se 
muestra en la tabla IV. Los números de las columnas segunda y ter- 
cera de la tabla II pueden ajustarse casi exactamente a las diferen- 
cias de un error único E, que se muestra en la tabla IV, si tomamos 
E = —54, y entonces la tabla IV queda 


0 

0 
0 0 

0 — 54 
0 — 54 

54 162 

— 54 108 

54 — 1162 
0 — 54 

0 54 
0 0 0 
0 


Así que a partir de solo las diferencias de la tabla II tenemos una 
fuerte evidencia de que su entrada 0,870 239 en x = 0,515 es un va- 


lor erróneo de la forma cos (0,515) + E con E = —54 y por consiguien- 
te que 
cos (0,515) = 0,870 239 + 54 unidades en la última cifra decimal 


0,870 293 


Este es en efecto el valor correcto, como se da en la tabla IL 
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Tabla V 


Ejercicio 1 


Localizar y corregir el error en la siguiente tabla de raíces cúbicas 


x Yx 


197 5, 818 648 

- 198 5,828 477 
_ 19 5, 838 272 
200 5, 848 035 
201 5,857 776 
202 5,867 464 
203 5,877 131 
204 5,886 765 
205 5,896 369 
206 5,905 941 


4.4.2 Errores inherentes 


Aun si no hay equivocaciones en las mediciones o en los cálculos cu- 
yos resultados se anoten en una tabla, los números tabulados tienen 
tendencia a estar erróneos, tanto por la exactitud limitada del proceso 
que los produce, como por el hecho de que se redondean con un núme- 
ro determinado de cifras decimales. En la Unidad 2, Errores y exacti- 
tud usamos el nombre de errores inherentes para todos los errores de 
este tipo. El propósito de esta sección es mostrar al lector cómo tales 
errores limitan la exactitud de los cálculos basados en una tabla de 
diferencias. Consideraremos principalmente los errores de redondeo, 
pero las mismas consideraciones se aplican a todos los errores inhe- 
rentes. Como ejemplo para ilustrar el efecto de los errores de redondeo 
sobre una tabla de diferencias, consideramos las diferencias de una 


2 
tabla de dos cifras de valores de 5 Para x = 31, 32,..., 37. 


" Diferencias 
dá 3 
Primera Segunda Tercera Cuarta 

31 320,33 
2 305 e WE 
33 363,00 . 0,66 y + 0,02 

22,33 + 0,01 —0,03 
34 385,33 23.00 0,67 000 7 0,01 0.00 
35 408,33 2367 0,67 001 —0,01 ? 
36 432,00 2433 0,66 , 
37 456,33 ? 


e e ' 
Como En es un cuadrático, esperamos que las segundas diferencias 


sean constantes y las diferencias de orden superior sean cero; pero 
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Ejercicio 1 
(5 minutos) 


4.4.2 


Discusión 
* 


Tabla 1 


por razón de los errores de redondeo los valores tabulares no son los 


x , : z 
valores exactos de + y en consecuencia las diferencias de los valores 


tabulares no son exactamente constantes en la columna de las segun- 
das diferencias y no son exactamente cero en las últimas columnas 
de diferencias. Cuando se usan tablas de diferencias, es importante 
saber si las diferencias de la tabla reflejan una propiedad real de la 
función tabulada como las segundas diferencias 0,66 y 0,67 de la ta- 
bla 1, o si son solamente diferencias de error de redondeo, como las di- 
ferencias tercera y cuarta de la tabla 1 En la tabla lI la distinción es 
obvia, ¿pero, es tan obvia en la tabla siguiente? 


Diferencias 


x Fx) 


Primera Segunda Tercera Cuarta Quinta 


0,25 0,247 404 


0,26 0257081 5 -27 . 

0,27 0266731 g0e 25 Z - ñl 
0,28 0,276356 ¿39 29 ' E 
0,29 0,285952 gero  —28 A. 2 : 
0,30 0,295520 ge 2 7. -l ; 
0,31 0,305059 g30s 31 " 1 - 
0,32 0314567 gave 32 , IS 
0,33 0,324043 yy 32 ha 2 

0,34 0,333487 gq) 33 


0,35 0,342 898 


Las diferencias decrecen en magnitud al principio, luego se incremen- 
tan a medida que el orden de las diferencias se hace mayor. ¿Puede 
decir si estos cambios representan una tendencia sistemática, o si se 
debe puramente a errores de redondeo? Tales preguntas nos conducen 
cerca de la frontera entre nuestro tema presente, diferencias finitas, 
y el problema de desenmarañar tales efectos sistemáticos a partir de 
los efectos de fluctuaciones erráticas tales como errores de redondeo. 
El enfoque estadístico de este tipo de problema se puede tratar en las 
unidades sobre estadística, pero aquí son suficientes los métodos 
de la Unidad 2, Errores y exactitud. Usaremos estos métodos para cal- 
cular límites de los errores en la tabla de diferencias producidos por 
errores de redondeo en los valores tabulados. Esto nos conducirá a 
deducir que una entrada cualquiera de la tabla de diferencias que sea 
mayor que este límite de error puede aparecer por lo menos en parte 
por otra causa diferente que por error de redondeo. Esto nos ayudará 
a decidir sobre el grado del polinomio con el cual podemos aproximar 
razonablemente la función tabulada. 


Para el límite de error absoluto en A»f(xx*) = flxx+1) — f(x») agre- 
gamos los límites de error absoluto de redondeo en los términos sepa- 
rados fíxx+1) y f(x»), de acuerdo con las reglas descritas en la Uni- 
dad 2, Errores y exactitud. El límite de error absoluto de redondeo en 
flxx+1) — flor) es por consiguiente ¿3 + 3 = 1 unidad en la última 


cifra decimal. Para el límite de error absoluto de redondeo en 
Ar. (xi) = Anflxi+ 1) — Apt xi) 


* La ilustración se toma de P. Henrici, Elements of Numerical Analysis, (John Wiley 1964). 


55 


MB4.4.2 


Tabla Il- 


(continúa en la página 36) 


MB4.4.1/4.4.2 
Solución 4.4.1.1 Solución 4.4.1.1 


La tabla de diferencias de los valores dados es 


5,818 648 

9829 
5,828 477 - 34 

9795 +2 
5,838 272 —32 

9763 +10 
5,848 035 

9741 —31 
5,857 776 —53 

9688 4-32 


5,867 464 


9667 
5,877 131 

9634 +3 
5,886 765 -30 

9604 3 
5,896 369 y 

9572 
5,905 941 


Aquí las terceras diferencias de la función tabulada no son constan- 
tes, pero el esquema es consistente con un error de E = +10 que con- 
tribuiría a la tabla de diferencias así: 


+10 
+10 
+10 —30 
+10 —20 
—10 +30 
+10 
—10 


Así que la tabla de diferencias indica que el valor 5,857 776 del ápice 
del triángulo tiene un error de +10 y que el valor correcto es 5,857 766. 
ES] 


(viene de la página 55) 


de nuevo agregamos los límites de error absolutos en los términos se- 
parados Arfílxi+1) y Anfílxr). El cálculo que acabamos de hacer 
muestra el límite del error absoluto de redondeo de A,f(x»*) que 
es 1 unidad en la última cifra decimal y que es el mismo en el de 
Arflxx+1) por un cálculo análogo. Por la ley de adición para lími- 
tes de error absoluto, el límite de error absoluto de redondeo en 
Ajf(xx) es por tanto 1 + 1 = 2. De la misma manera el límite de 
error absoluto de redondeo de Ajf(xs) es 2 + 2 = 4, el de Af f(x») 
es 4 + 4 =8 y en general, 


el límite de error absoluto de redondeo de A;f(x*) es 2”7' 


unidades en la última cifra decimal. 
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Ejercicio 1 

Calcular los errores producidos en una tabla de diferencias por redon- 
deo, si los errores de redondeo en los valores tabulares son alternati.- 
vamente +3 unidades y —¿ unidades en la última cifra decimal y 


comparar con el límite de error absoluto de A;f(xx), dado antes. 
[Es] 


Volviendo al ejemplo considerado en la tabla II podemos ahora poner 
a prueba el significado de la tendencia de los números de las colum- 
nas de diferencias. Por ejemplo, el límite de error absoluto de redon- 
deo de las quintas diferencias es 2% = 16; en consecuencia, cuales- 
quiera sean las quintas diferencias cuyo valor numérico sea menor 
que 16 puede deberse puramente a error de redondeo. Concluimos que 
las quintas diferencias tabuladas no dan una evidencia firme acerca 
de la naturaleza de la función tabulada. Buscamos los de las columnas 
de diferencias en las cuales las entradas puedan no ser atribuidas en- 
teramente a errores de redondeo. 


Ejercicio 2 


¿Cuáles de las diferencias de la tabla II son apreciablemente diferen- 
tes de cero (i.e. no pueden atribuirse enteramente a redondeo) ? E 


De la respuesta del último ejercicio podemos sacar la conclusión de 
que no tiene caso ajustar los valores tabulados de la función de la ta- 
bla II con un polinomio de grado superior al del cuadrático. Si tratá- 
ramos de usar algún grado más alto deberíamos estar ajustando el 
polinomio a los errores de redondeo y no a la función misma. Esto 
implica que para interpolar en la tabla II, una fórmula de interpolación 
cuadrática, tal como la de Gregory-Newton cuadrática, dada por la 
ecuación (4) de la Sección 4.3.2 da la más alta exactitud que esta ta- 
bla puede dar. 


Ejercicio 3 


¿Cuál es el grado de interpolación polinómica óptima para la interpo- 
lación de la tabla siguiente? Uselo para hallar tan (1,4975). 


Xx tan x 


1,495  13,1680 
1,496  13,3447 
1,497  13,5262 
1,498 13,7127 
1,499  13,9043 
1,500  14,1014 
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Ejercicio 1 
(3 minutos) 


Discusión 
* 


Ejercicio 2 
(2 minutos) 


Discusión 
* 


Ejercicio 3 
(3 minutos) 


Solución 1 


Primeras Segundas Terceras Cuartas 


dif. dif. dif. dif 
1 
z 
+) 
1 
> 2 
2 = 
Í . -2 E 8 
3 E 2 : -3 
2 1 —á 
j 3 
ad E 
2 


Se ve que las diferencias de los errores se doblan, en este caso, en ca- 
da vez que se hace la diferencia. Vemos entonces que con esta esco- 
gencia de errores de redondeo obtenemos el límite de error absoluto de 
redondeo de Ar fíx*), i.e. 271. A 


Solución 2 


La primera y la segunda diferencias tienen significado definido y son 
diferentes de cero. El límite de error absoluto de redondeo de las ter- 
ceras diferencias es 2? = 4 así que cualesquiera sean las terceras 
diferencias cuyo valor numérico sea menor que 4, esto puede deberse 
a error de redondeo. Así que deberíamos considerar las terceras dife- 
rencias como no apreciablemente diferentes de cero. a 


Solución 3 


Diferencias 
Xx tan x 
Primera Segunda Tercera Cuarta 


1,495  13,1680 


1,496  13,3447 e 48 s 

1,497 13,5262 ¡26 s0 $ 
1,498 18,127 jo 3 p 

1,499 — 13,9043 ¡971 55 


1,500  14,1014 


Para las segundas diferencias el límite de error absoluto es 2 unidades 
en la última cifra decimal; por consiguiente, las segundas diferencias 
tabuladas son significativas. Para las terceras diferencias el límite de 
error absoluto es 4; por consiguiente las terceras diferencias tabula- 
das no son significativamente diferentes de cero. Entonces la fórmula 
de interpolación que da la exactitud más aceptable sin demasiado tra- 
bajo, es la que se obtiene tratando las terceras diferencias como cero 
y es por tanto, de segundo grado. 


Por la fórmula de Gregory-Newton, 


0(0 — 1 
tan (1,4975) = 13,5262 + 0,18650 + 0,0051 ( > ont = $ 
= 13,5262 + 0,0932 — 0,0006 
= 13,6188 Ll 
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Solución 1 


Solución 2 


Solución 3 


4.43 Resumen 


En esta última parte hemos mostrado cómo las equivocaciones llevan 
a una tabla de diferencias de la forma 


E 

E 
+E —3E 

E —2E 
—E +3E 

E 
Y 


que con frecuencia facilita la detección y corrección de tales equivo- 
caciones. 


Por otra parte, los errores de redondeo de +3 en la última cifra deci- 
mal llevan a límites de error para diferencias sucesivas de 1, 2. 4.8, 
respectivamente. Por tanto, si hallamos las magnitudes de las dife- 
rencias más allá del orden r sin exceder a 2”-! en la última cifra de- 
cimal, podemos concluir que la función es un polinomio de grado r, con 
errores de redondeo, o, para precisar más, que no tiene caso tratar de 
representar la función por un polinomio de grado superior a r. Los mé- 
todos de diferencias finitas llevan a métodos de interpolación y extra- 
polación de cualquier función continua y de deducción de fórmulas 
generales. 


Esto concluye nuestra consideración de diferencias finitas, un tema 
de vital importancia en el cálculo práctico, que también arroja una luz 
interesante sobre el comportamiento de las funciones en general e in- 
troduce ideas que se usarán en el trabajo futuro sobre sucesiones, in- 
tegración y diferenciación. 
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Resumen 
* «$ 


Reconocimientos 


Expresamos nuestra gratitud a las siguientes fuentes por el material 
usado en este texto por correspondencia: 


Textos 


John Wiley and Sons Ltd, Handbook of Engineering Fundamentals 
por O. W. Eshbach, 1966, y Elements of Numerical Analysis por P. 
Henrici, 1964; 

Cambridge University Press, Science and Music por Sir James Jeans, 
1961; 

The Rt. Hon. P. J. Noel-Baker, The Arms Race publ. por Oceana Pu- 
blications 1956. 


Ilustraciones 


The Mansell Collection, fotografías de Karl W. T. Weierstrass y Jo- 
seph-Louis Lagrange. 
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M 100 — UNIDADES DEL CURSO BASICO DE MATEMATICAS 


00-250 0O0*>0NA 


Funciones 

Errores y exactitud 
Operaciones y morfismos 
Diferencias finitas 

NO HAY TEXTO 
Desiguales 

Sucesiones y límites 1 
Computación I 

Integración I 

NO HAY TEXTO 

Lógica 1 — Algebra de Boole 
Diferenciación Il 
Integración II 

Sucesiones y límites II 
Diferenciación II 
Probabilidad y estadística I 
Lógica 11 — Prueba 
Probabilidad y estadística II 
Relaciones 

Computación II 
Probabilidad y estadística II 
Algebra lineal I 

Algebra lineal II 
Ecuaciones diferenciales 1 
NO HAY TEXTO 

Algebra lineal HI 

Números complejos I 
Algebra lineal IV 

Números complejos II 
Grupos I 

Ecuaciones diferenciales II 
NO HAY TEXTO 

Grupos II 

Sistemas numéricos 
Topología 

Estructuras matemáticas 
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